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The  first  portion  of  this  document  deals  with  the 
synthesis  of  spherical  mechanisms  and  is  given  with  the 
formulation  of  the  necessary  tools.   These  tools  are  the 
curvature  transformation  and  the  inversion  principle  which 
allow  closed-form  design  of  multi-link  mechanisms  for 
three,  four,  and  five  multiply  separated  positions. 
Twelve  solvable  problems  are  tabulated  with  an  emphasis  . 
on  the  reduction  of  parameters  for  motion  specification. 
The  curvature  transformation  is  extended  to  include  geo- 
metric or  time  base  motion  problems.   More  than  20  detailed 
design  procedures  are  presented  for  the  spherical  four-, 
five-,  and  six-link  mechanisms. 

The  second  portion  of  this  document  deals  with 
the  analysis  of  spherical  mechanisms  and  is  presented  in 


vi 


terms  of  two  complementary  formulations.   The  first  pro- 
vides for  trans! ational  and  angular  evaluation  of  multi- 
link  mechanisms.   This  analysis  is  established  by  means 
of  two-link  chains  known  as  dyads.   The  second  treats  the 
input-output  equation  algebraically.   This  format  yields 
the  evaluation  of  multiple  closure  branches  simultaneously. 

All  of  the  synthesis  and  analysis  tools  are  shown 
to  be  reducible  for  use  in  planar  mechanism  design. 


VI  1 


.-♦^v.-' 


'■  (* 


CHAPTER  I 
INTRODUCTION 


The  field  of   kinematics  is  becoming  a  major  con- 
tributor toward  advances  in  industrial  automation.   Ma- 
chines are  required  to  simulate  motions  of  human  operations 
as  well  as  motions  beyond  human  capability.   It  is  demanded 
that  these  functions  be  accurate,  reliable,  and  efficient. 
Such  demands  have  created  growing  interest  in  the  field 
of  machine  design.   Special  spatial  motion  problems  such 
as  transfer  mechanisms  in  automatic  processing  equipment, 
function  generation  between  nonparallel  axes,  and  path 
generation  on  a  sphere  cannot  be  solved  by  planar  mecha- 
nisms.  The  relative  importance  of  the  spherical  system 
can  be  obtained  by  considering  that  perhaps  75%  of  the 
industrial  linkage  problems  can  be  satisfied  by  planar 
mechanisms,  15%  by  spherical  mechanisms,  and  the  remainder 
by  more  complex  spatial  mechanisms. 

Coplanar  synthesis,  or  design  of  planar  linkages, 
has  reached  a  high  level  of  sophistication,  usefulness, 
and  completeness— a  development  which  has  occurred  chiefly 
over  the  past  15  years.   Analytical,  graphical  and  special 
design  procedures  are  now  widely  available  to  create  four- 

1 


five-  and  six-bar  mechanisms.   Proper  application  of  the 
four  tools  of  synthesis  (inversion,  angular  cognates,  path 
cognates,  and  curvature  theory)  allows  the  design  of  the  geared 
five-bar,  Stephenson  six-bar,  and  Watt  six-bar  for  a  wide 
range  of  motion  problems  basic  to  industrial  application. 
Since  the  spherical  linkage  has  properties  which 
are  distinct  from  those  of  coplanar  linkages  and  repre- 
sents perhaps  15%  of  potential  industrial  application,  in- 
creased competence  by  the  designer  is  warranted.   Fortu- 
nately, much  useful  work  on  the  spherical  four-bar  has  oc- 
curred, although  its  general  level  is  approximately  that 
of  coplanar  synthesis  15  years  ago.   In  particular,  few 
design  procedures  have  been  evolved  for  linkages  more  com- 
plex than  the  spherical  four-bar,  such  as  the  spherical 
bevel  geared  five-bar,  and  equivalent  six-bar  linkages. 
The  development  of  the  spherical  tools  of  synthesis  would 
fulfill  this  need  for  the  spherical  system.   A  unified 
theory  of  spherical  kinematics  would  not  limit  itself  to 
the  suggested  15%  application.   Rather  it  would  provide 
common  tools  for  both  the  plane  and  the  sphere.   Since 
planar  geometry  is  a  proper  subset  of  spherical  geometry, 
the  concepts  would  cover  90%  of  application  problems.   Such 
a  unified  theory  may  prove  to  be  the  key  to  a  total  design 
package  with  the  aid  of  dual  numbers.   Future  efforts  may 
extrapolate  spherical  kinematics  to  provide  unified  ana- 
lytics of  planar,  spherical  and  spatial  motions. 
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Hence,  the  objective  of  this  document  is  to  de- 
velop general  techniques  which  would  make  the  spherical 
system  equally  as  accessible  to  the  designer  as  the  co- 
planar  system.  Every  effort  has  been  expended  to  insure 
that  the  results  are  immediately  useful  and  valuable  to 
the  designer. 

Generally  the  approach  to  spherical  design  has 
been  to  extend  well-established  methods  for  synthesis  in 
the  plane  to  the  sphere.   For  instance,  J.  Denavit  and 
R.  S.  Hartenberg  [1]  generalized  F.  Freudenstei n ' s  approx- 
imate synthesis  of  planar  four-bar  linkages  to  generate  . 
functions  by  means  of  the  spherical  four-bar.   Also,  J.  T, 
Wilson  [2],  H.  L.  Johnson  [3],  and  J.  R.  Zimmerman  [4]  ob- 
tained design  equations  for  function  generation  using  a 
rotation  matrix. 

..  C.  H.  Suh  and  C.  W.  Radcliffe  [5]  specified  mul- 
tiple positions  on  a  path  in  the  plane  by  a  displacement 
matrix.   Following  this  they  extended  the  procedure  to 
the  sphere  using  a  rotation  matrix.   These  contributions 
emphasized  the  importance  of  matrix  methods  for  synthesis 
introduced  by  J.  Denavit  and  R.  S.  Hartenberg. 

General  theories  for  finitely  and  inf ini tesimal ly 
separated  position  motion  of  a  rigid  body  in  space  have 
been  developed  by  B.  Roth  [6,  7,  8,  9,  10]  and  G.  N. 
Sandor  [11].   These  are  major  contributions  in  the  field 
of  synthesis  of  spatial  motions.   In  [9,  10],  P.  Chen  and 
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B.  Roth  developed  a  general  theory  which  enables  the  synthe- 
sis of  mixed  displacements  (i.e.,  finite  and  infinitesi- 
mal combined)  of  a  rigid  body  in  space.   In  [11]  G.  N. 
Sandor  developed  a  general  method  for  the  synthesis  of 
spatial  mechanisms  using  quaternions. 

The  theory  of  coplanar  synthesis  has  been  studied 
exhaustively.   The  study  and  development  of  a  correspond- 
ing theory  for  synthesis  of  spherical  motion  has  been  much 
more  limited.   G.  Dittrich  [12],  H.  J.  Kamphuis  [13],  and 
K.  E.  Bisshopp  [14]  have  studied  instantaneous  spherical 
kinematics. 

The  design  of  the  planar  geared  five-bar  for  . 
function  generation  was  treated  by  Oleksa  and  Tesar  [15]. 
The  extension  of  this  philosophy  to  the  spherical  five-bar 
was  this  author's  first  introduction  into  spherical  kine- 
matics.  The  foundation  of  this  extension  was  based  upon  the 
unified  theory  of  spherical  trigonometry  by  Duffy  and 
Rooney  [16].   This  theory  was  generated  from  an  understand- 
ing of  the  work  by  Todhunter  and  Leatham[17].   Dowler, 
Duffy,  and  Tesar  [18-20]  combined  the  theory  of  spherical 
trigonometry  and  the  planar  Burmester  theory  to  yield  the  ■ 
spherical  curvature  transformation.   This  tool  gave  the 
author  the  capability  to  extend  the  function  generation 
problem  to  the  sphere  [21,22].   Duffy,  Riddle,  and  Tesar  [23] 
then  established  a  displacement  analysis  technique  for  the 
spherical  five-bar. 


Eschenback  and  Tesar  [24]  expanded  the  synthesis 
concept  of  time  dependent  motion  for  the  planar  system. 
Myklebust  and  Tesar  [25,26]  extended  this  effort  to  yield 
time  dependent  motion  coefficients  for  the  curvature  trans- 
formation to  establish  a  uniform  approach  to  treat  this  new 
problem.   A  generalized  computer  program,  MECSYN,  was  de- 
veloped for  the  synthesis  of  multi-link,  planar  mechanisms 
with  time  state  motion  specification.   A  collection  of  de- 
sign procedures  [27]  by  Tesar  for  finitely  separated  motion 
specification  has  been  useful  for  academic  instruction.   A 
major  portion  of  the  included  work  will  be  devoted  to  estab- 
lishing similar  design  procedures  for  the  spherical  system. 

Analysis  provides  an  evaluation  tool  of  the  use- 
fulness of  a  synthesized  mechanism.   Gilmartin  and  Duffy 
[28]  give  a  means  to  identify  the  type  of  a  four-bar 
mechanism.   Freudenstein  and  Primrose  [20]  have  studied 
the  motion  of  spherical  and  spatial  five-bars  for  closure 
conditions.   The  most  useful  analysis  tool  for  planar 
mechanisms  has  been  established  by  Pollock  and  Tesar  [30] 
for  dyad-based  machinery.   This  philosophy  will  be  con- 
sidered for  the  analysis  of  spherical  motion  in  the  second 
portion  of  this  work. 

It  will  be  the  object  of  this  work  to  present 
spherical  synthesis  and  analysis  tolls  in  a  form  usable 
to  industrial  designers. 


CHAPTER  II 
SYNTHESIS  CONCEPTS 


The  creative  design  of  machine  systems  needs  to 
be  based  upon  a  solid  foundation  of  analytical  understand- 
ing.  The  technological  needs  of  today  have  progressed 
beyond  the  trial  and  error  evolution  of  past  machine  design 
practice.   The  development  of  sophisticated  analytical 
formul ations, coupl ed  with  the  practical  convenience  of  com- 
puters, has  given  the  designers  tools  to  accomplish  their 
task  with  greater  efficiency.   The  synthesis  concepts  in- 
cluded within  are  presented  with  enough  detail  so  that 
they  might  be  useful  to  all  engineers  interested  in  de- 
signing machines  with  spherical  constraints. 

It  is  essential  that  the  designer  recognize  the 
number  of  parameters  {S}  that  must  be  specified  to  satisfy 
the  problem  requirements.   If  the  designer  selects  a  mech- 
anism with  too  few  available  parameters,  {q},   he  will  be 
severely  limited  in  the  quality  of  results.   On  the  other 
hand,  if  the  designer  selects  a  mechanism  with  too  many 
available  parameters,  the  range  of  solutions  will  be  too 
vast  for  a  meaningful  evaluation  of  possible  solutions. 


A  basic  criterion  of  the  design  process  should  be  that  the 
available  parameters  {q}  exceed  or  equal  the  specified 
parameters  {S}  or  - 


{q}  >   {S} 


(II-l) 


The  synthesis  problem  involves  the  satisfaction  of 
a  number  {F}  of  conditional  functions  in  a  limited  number 
{i  +   1)  of  system  positions.   Suppose  the  system  is  to 
satisfy  the  condition  f(x,y,z)  =  0  by  the  approximating 
function  g(x,y,z)  =  0  as  in  Fig.  II-l.   The  function 
f(x,y,z)  may  be  a  general  spatial  curve  or  it  may  be  spe- 
cialized to  spherical  motion  by  the  constraint 


2    2    2 

X  +  y  +  z   =  constant 


(II-2) 


The  function  may  be  further  specialized  to  planar  motion 
by  restricting  one  of  the  parameters  (x,y,z)  to  be  a  con- 
stant. 

Closed  form  synthesis  techniques  allow  a  maximum 
of  five  design  positions.   This  does  not  restrict  the  ac- 
tual number  of  intersections  of  the  two  curves,  but  rather 
the  number  of  positions  free  to  the  desi gner's  specif ica- 
tion.   These  positions  may  be  either  finitely  or  infin- 
itesimally  separated.   Inf i ni tesimal ly  separated  positions 
yield  the  ability  to  match  slope,  curvature,  etc.  up  to 
the  fourth  order.   After  the  choice  of  design  positions, 
the  designer  must  consider  whether  it  is  necessary  for 
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f(x,y,2) 


g(x»y,2) 


Figure    II-l 


g(x,y,z)  to  be  spatial,  spherical  or  planar  to  meet  the 
specified  design  positions.   It  should  be  noted  that 
spatial  mechanisms  include  spherical  and  planar  mechanisms 
as  a  proper  subset.*   Planar  mechanisms  are  also  a  subset 
of  spherical  mechanisms.   With  this  in  mind,  the  designer 
could  always  choose  spatial  specifications  to  satisfy  the 
problem  function  f(x,y,z).   Spatial  mechanisms  have  many 
more  parameters  to  be  considered  and  are  often  unneces- 
sarily complex.   In  addition,  there  are  no  well-established 
synthesis  techniques  available  for  general  spatial  mech- 
anisms.  References  include  spatial  analysis  routines  that 
may  be  considered  for  trial  and  error  design,  but  this  is 
recommended  only  as  a  last  resort.   If  f(x,y,z)  is  a 
general  spatial  function,  it  would  be  most  desirable  to 
approximate  it  with  a  spherical  function  g(x,y,z).   Spher- 
ical mechanisms  add  little  complexity  to  the  analytics, 
yet  allow  a  type  of  spatial  motion.   It  will  also  be  seen 
that  spherical  and  planar  mechanisms  may  share  the  same 
body  of  analytics  and  add  no  extra  effort  in  use.   This 
work  will  continue  to  emphasize  spherical  synthesis  but 
the  same  philosophy  and  reduceabl  e analytics  follow  directly 


♦Spatial  mechanisms  have  no  restrictions  on  the 
rotation  axes.   Spherical  mechanisms  have  all  rotation  axes 
intersecting  at  common  point,  the  origin  of  the  sphere. 
Planar  mechanisms  have  all  rotation  axes  parallel,  or  in- 
tersecting at  infinity. 
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to  the  plane.   We  can  choose  only  a  limited  number  of  de- 
sign positions  (x,  y. ,  z^).   ^=  0,1,2,  ...   to  repre- 
sent f(x,y,z).   Should  more  than  one  function  (F  =  1,2,3, 
...)   need  to  be  satisfied  simultaneously,  the  total  number 
of  specified  parameters  would  be 


S  =  F£ 


(II-3) 


Note  that  this  result  does  not  include  the  1=0    position. 
Mechanism  constraints  can  satisfy  di spl acement,  but  normally 
not  absol ute,  val ues .   This  in  no  way  weakens  the  synthesis 
approximation,  since  any  functional  reference  may  be  chosen, 
Ideally,  the  synthesis  process  will  not  utilize  all 
the  system  parameters  {q}  to  satisfy  the  specifications  {S}, 
The  remaining  parameters 


F£ 


(II-4) 


can  be  used  to  meet  optimization  criteria  not  included  in 
the  specifications.   A  synthesized  mechanism  will  have 
zero  error  at  the  Z   +    1  design  positions,  but  no  other 
characteristics  are  guaranteed.   There  are  over  a  dozen 
conditions  to  be  met  after  synthesis.   All  of  these  can 
be  evaluated  through  present  analysis  procedures.   The 
designer  could  be  left  with  thousands  of  these  evalu- 
ations, a  number  much  too  great  to  treat  realistically. 
Optimization  is  the  answer,  and  it  is  now  coming  to  the  • 
surface  for  planar  mechanisms.   Concise  and  efficient 
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synthesis  and  analysis  procedures  are  the  necessary  keys 
to  this  type  of  undertaking.   The  included  work  should 
make  a  major  contribution  in  this  arena  for  spherical 
mechanisms.   The  number  of  optimization  parameters  {p} 
should  be  in  the  range  1  1  P  1  3  for  good  design. 

Constraint  Systems 


The  designer  has  a  wide  choice  of  mechanisms  to 
meet  the  requirements  of  a  given  problem.   The  control 
system  contains  a  few  parameters,  usually  one  or  more  of 
which  is  completely  adjustable.   The  cam  system  contains 
an  infinite  collection  of  nonadjustabl e  parameters.   A 
linkage  is  neither  adjustabl e,  nor  does  it  contain  an  in- 
finite collection  of  parameters.   The  designer  must  know 
whether  the  linkage  contains  enough  parameters  {q}  to 
satisfy  his  problem.   This  work  will  build  a  family  of 
solutions  based  upon  the  four-bar,  the  geared  five-bar, 
and  the  six-bar  mechanisms.   As  in  the  plane,  it  is  felt 
that  these  are  the  building  blocks  of  linkage  design  and 
a  very  large  percentage  of  solutions  can  be  derived  from 
their  motions. 

The  spherical  four-bar  mechanism  (Fig.  II-2)  is 
the  most  elementary  device  of  this  study  and  contains  the 
fewest  system  parameters  {q}.   If  we  were  to  ask  the  ques- 
tion, "What  is  the  fewest  number  of  parameters  necessary 
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Figure    1 1-2 
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to  specify  the  geometric  position  of  this  device?",  the 
answer  could  be  given  in  any  number  of  terms.   The  most 
immediate  answer  is  derived  by  considering  the  location 
of  each  pin  joint  A,B,C,D.   If  the  condition  of  a  spheri 
cal  mechanism  is  satisfied,  then  the  axes  of  each  pin 
joint  passes  through  the  origin  0.   The  arc  between  any 
two  pin  joints  is  a  great  circle  segment.   If  we  arbi- 
trarily select  a  unit  sphere  for  convenience,  such  that 


z2  =  1  _(x2  +  y2j, 


(II-5) 


only  two  parameters  would  be  necessary  to  describe  the 
position  of  each  pin  joint.    A  total  of  q  =  8  parameters 
is   necessary  to  complete  this  description.   The  special 
case  of  a  prismatic  slider  can  be  simulated  on  the  sphere 
by  a  90°  crank.   If,  in  Fi g.  II-2, arc  CD  is  90°,  then  point 
C  will  move  on  a  great  circle  path.   This  is  the  equiva- 
lent of  linear  motion  on  the  plane. 

The  spherical  five-bar  mechanism  (Fig.  II-3a)  has 
two  degrees  of  freedom.   Any  two  links  of  the  chain  can 
be  driven  independently  to  constrain  the  motion.   Alter- 
natively, gear  trains  can  be  introduced  between  two  non- 
adjacent  links  to  yield  a  one  degree  of  freedom  device. 
Motions  derived  from  the  five-bar  mechanism  family  are 
inherently  more  complex  than  those  of  the  four-bar  mech- 
anisms.  It  is  expected  that  a  geared  five-bar  mechanism 
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(b)  FIVE-BAR     I     (fix    link  I) 

(c)  FIVE -BAR    2    (fix    link  2) 

(d)  FIVE-BAR    3    (fix    link  3) 


Figure    1 1 -3 
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will  solve  some  problems  that  a  four-bar  would  not.   Note 
that  there  are  five  pin  joints  in  the  five-bar,  so  that 
10  parameters  are  necessary  to  describe  its  location.   In 
order  to  predict  the  position  of  all  links  relative  to  a 
new  input  position,  two  additional  parameters  must  be  con- 
sidered.  The  gear-train  values  M  and  N  (Fig.  II-3b,c,d) 
are  necessary  to  describe  the  input  dependence  of  the 
coupled  1  inks.   Therefore,  a  total  of  q  =  12  parameters  are 
necessary  for  the  geared  five-bar  mechanisms. 

The  Watt  six-bar  mechanism   (Fig.  II-4)  has  seven 
pin  joints  so  that  q  =  14  parameters  is  necessary  to  de- 
scribe its  position.   The  Watt  1  and  Watt  2  mechanisms  are 
distinct  inversions  of  the  general  Watt  chain  that  are 
found  by  fixing  links  1  and  2  respectively.   The  reader 
should  recognize  that  the  motion  in  link  5,  for  example,.  . 
is  more  complex  in  the  Watt  1  than  that  in  the  Watt  2.   It 
is  true  that  the  motions  of  six-bars  with  two  fixed  pivots 
are  more  complex  than  those  with  three  fixed  pivots. 

The  Stephenson  six-bar  mechanism  (Fig.  II-5)  also 
requires  q  =  14  parameters.   It  has  three  distinct  in- 
versions found  by  fixing  links  1,  2,  and  3. 

Solvable  Problems 
With  a  view  of  the  components  in     . 


p  =  q  -  Fi 


(II-6) 
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Figure    II-4 
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Figure    II-5 
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we  are  now  ready  to  consider  the  range  of  solvable  problems 
and  the  selection  of  the  proper  mechanism  for  each  type. 
This  range  of  problems  can  be  most  efficiently  presented 
by  the  use  of  tables  with  consistent  nomenclature. 

Let 


Z   =  a, 1,2, 

q 

F 

a,b,c ,  ... 

m  =  T,2,3, 
m  =  1  ,2,3. 

e  : 

E{u,v) 

2:(u.v) 


^m 


5i 


the  multiply  separated  position 
counter 

the  number  of  parameters  left  for 
optimization  (p,  =  q  -  F, ) 

total  available  parameters  in  the 
constraint  system 

number  of  functions  to  be  specified 
simul taneously 

pin  connectors  in  the  solvable 
probl em  statement 

pin  connectors  in  the  linkage 
system 

counter  for  links  in  the  solvable 
problem  statement 

counter  for  links  in  the  linkage 
system 

point  in  the  problem  tracing  a  de- 
sired point  path  with  curvilinear 

coordinates  (U  ,V  ) 
e  e' 

moving  shel 1 

reference  shel 1 

absolute  angular  motion  parameter 
in  the  solvable  problem  statement 

relative  angular  motion  parameter 
in  the  solvable  problem  statement 
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[m, ,  m^,  m  ]         -  available  system  links  to  repre- 
'   ^        3  sent  links  in  the  problem  state- 

ment 

p  -  pin  joint  available  to  represent 

a  general  path  point  e 


The  tables  on  the  following  pages  are  intended  as  design  sheets 
for  each  of  the  specified  problems  to  assist  the  designer 
in  the  first  stages  of  the  design  process.   The  tables  are 
oriented  on  the  base  of  finitely  separated  position  (FSP) 
problems  but  are  extrapolated  to  include  guides  for  the 
use  of  multiply  separated  positions  (MSP)  and  time-state 
synthesis  (TSS).   The  distinctions  of  FSP,  MSP,  and  TSS  . 
will  be  made  clear  in  the  curvature  transformation  dis- 
cussion. 

Synthesis  Tools 

..;_■    Planar  kinematics  have  four  well-established  tools 
of  synthesis. 

1 .  path  cognates 

2.  angular  cognates 
3  .   i  n  V  e  r  s  i  0  n 

4.   curvature  transformation 

These  are  fully  discussed  in  the  literature  [27]  for  planar 
mechani  sms  but  have  major  distinctions  for  spherical  synthesis 
Both  the  path  and  angular  cognates, that  are  so  helpful  in 
planar  design, are  based  upon  properties  of  parallelograms. 
Unfortunately,  there  are  no  parallel  great  circles  on  the 
sphere,  and  these  two  tools  must  be  excluded  from  the  begin- 
ning. 
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Spherical  Chain 

To  appreciate  the  effect  of  colaminar  inversion 
we  can  first  study  the  spherical  three-link  chain  (Fig. 
II-6).   The  motion  specification  set  {a.B,Y}^£  for  link 
m  in  position  I   with  reference  r  is  a  set  of  Euler  rota- 
tions about  the  X,y',z"   axes  (Fig.  1 1-7  ),  respecti  vely . 
From  the  sine,  sine-cosine  and  cosine  laws  for  the  spheri 
cal  triangle  abb',  we  write 

sin  Bp  ""  ^^'"  ^^  ^^""  ®l 


sin  a2  ^°^    ^2   ^  ^^""  ^^  ^°^  ^1 


(II-7) 


cos  ay    cos  ^2    ^    ^°^   ^^ 


Solving  for  B„  from  Eq.  (II-7), 


.  -1 


'2i 


dx 


-^  Csin"'  (sin  Tb  sin  9^)]^       (II-8) 


where  i   ^  MSP  counter 
k  ->-   ISP  counter 


T  ->■  any  geometric  or  time  parameter 


'■•V^     -p.'T'T.T 


33 


Figure    II-6 
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-y\-y" 


Figure    II-7 
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The  result  of  Eq.  (II-8)  for  k  =  0  will  yield  ^2  "i"  ^^^ 
first  or  fourth  quadrants.   This  is  sufficient  as  long  as 
ap  is  not  restricted  to  principal  quadrants.   From  Eq. 
(11-7),  ap  can  be  found  in  its  proper  quadrant  by  the  im- 
plicit expressions  for  sin  a^    and  cos  a^-      The  ratio  of 
the  last  two  expressions  in  Eq.  (II-7)  yields 


a 


Zi 


dy' 


tan 


-1 


sin  ab  cos  9 


1 


cos  a 


aF 


(II-9) 


The  two  parameters  a2j,  and  g^^  ^"^^  sufficient  to  describe 
the  location  and  higher  order  properties  of  point  b.   The 
third  Euler  angle  rotation  can  be  seen  (Fig.  II-6)  to  be 
the  sum 


Y 


2Z 


dy' 


(e2+p2) 


(11-10) 


The  angle  0^  is  assumed  to  be  given  as  a  function  of  e, 
with  its  appropriate  derivatives.   The  angle  pp  can  be 
represented  as  a  function  of  9-.  through  the  cosine  law  for 
thepolartriangleabb' 

0  =  sin  pp  cos  9-,  -  cos  p^  sin  9,  cos  ab 


(11-11) 
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or 


'2-1  "dx^ 


tan'   (tan  9 -j  cos  ab) 


(11-12) 


Eq.  (11-12)  can  be  differentiated  to  give  the  higher  order 
properties  of  Eq.  (11-10). 

Writing  the  sine,  sine-cosine,  and  cosine  laws  for 
spherical  triangle  ace*. 


sin  B3  =  X^  cos  0^  -  Y^  sin  6^ 


sin  ttj  cos  B3  =  -(X^  sin  9^  +  Y^  cos  9^) 


(11-13) 


cos  a^  cos  63  =  Z2 


where  X^  =  sin  be  sin  Q^^ 


Yj,  =  -(sin  ab  cos  be  +  cos  ab  sin  be  cos  d^^) 


Zp  =  cos  ab  cos  be  -  sin  ab  sin  be  cos  Qy 


Eq.  (11-13)  can  be  solved  to  find 
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CJT 


[sin"'  (X2  cos  e^  -  Y2  sin  e^)]j 


(11-14) 


a 


3£  -^^k 


tan 


_^  J  -X2  sin  e-i  -  Y2  cos  9 


This  can  be  extrapolated  to  a  general  n-link  chain  by  the 
following  expressions 


K.    -h     [^'"'^  <^n-l... 


dx 


32  ^°^  ^1  "  ^-1...32  ^^'"  ^l^^T 


(11-15) 


dx 


tan 


-1 


-X  ,   ,<,  sin  e,  -  Y^  ,  ^^   cos  9, 
n-1 . .  .32 1    n-l . .  .oc l_ 

^n-1...32 


The  final  rotation  angle  Y3  about  point  c  can  be 
expressed  as  the  sum  .  , 


Y3il=,^k 


[(93  +  03  +  p,)] 


3'-'T 


(11-16) 


as  depicted  in  Fig.  II-6.   For  the  spherical  right  triangle 
ace',  the  interior  angle  ^^   ^^"  ^^   expressed  as 
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tan  B- 


tan  d)-,  =  —. = 

*3   sin  a- 


(11-17) 


From  the  cosine  law  for  the  same  triangle, 


cos  ac  =  cos  a^  cos  6^ 


(11-18) 


Similar  to  Eq.  (11-12),  p_  can  be  expressed  implicitly  by 


tan  p_  =  tan  <})-  cos  ac 


(11-19) 


Substituting  Eqs.  (11-17,  18)  into  Eq.  (11-19),  we  find 


d        1 
^3i   "  ~~k  ttan"   (sin  6^  ^o*  '^S^^T 
dx 


(11-20) 


Consider  the  spherical  triangle  abc;  write  the  sine  and 
sine-cosine  laws 


X„  =  sin  ac  sin  a. 


Yp  =  -sin  ac  cos  a^ 


(11-21) 


where   X„  =  sin  ab  sin  6^ 


Y^  =  -(sin  be  cos  ab  +  cos  be  sin  ab  cos  9^) 
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Ratioing  Eq.  (11-21), 


3a  ~   ,    k 
dx 


/X, 


tan 


(11-22) 


Equation    (11-16)  can  be  extrapolated  for  a  n-link  chain 


Yn    =    ^      [9      +cr      +p]t 
'nil       .    k     '^    n  n        "^n   T 

dT 


(11-23) 


wh 


(k) 
ere  Q\    '    is  specified  and 


x\%       .  k 
dx 


tan 


-1/  ^23.  .  .n-1 
^23..  .n-1 


d  1 

Pf,^  =  — ^   [tan"   (sin  B^  cos  a„)]. 

dx 


'n'-"! 


The  computation  of  {a.B.Y}^^  an^l  {ct.B.Y}^^  for  the  three 
link  chain  has  been  validated  by  the  APL  function  LINK3  as 
listed  in  the  Appendix. 

Inversion 

The  motion  specification  set  {a.g.yJ^n  for  the 

III  X 

ttiree-lirik  chain  was  established  with  reference  to  link 
r.   We  will  see  in  the  design  procedures  that  it  is  de- 


sirable to  find  the  motion  specification  set  {a,3,Y}"„ 

'mil 
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with  respect  to  another  reference  link  n.   The  question 
might  be  posed,  "What  is  the  motion  specification  set 
{a,B»Y}oo  °'^  ^^^^   -^  viith    1  fixed  as  reference?"  To  answer 
this  question,  consider  the  transformation  between  the 
moving  reference  x,y,z  and  the  fixed  reference  X,Y,Z 
(Fig.  II-7) 


X 
Y 
Z 


=  M' 


m 


X 

y 

z 


(11-24) 


where  M   is  the  product  of  successive  rotation  matrices  of 
m        ^ 

Euler  angles  a  ,  3  ,  Y^  about  X,Y,Z  respectively,  i.e.. 


M' 


m 


(»ii):  (^2):  t^n': 


(^2,>: 


<^3i': 


'^22>: 


(>32>; 


(^23'I 


'^33'm 


(11-25) 


where 


(^ll^m  =  "^  ^m  "^  Y 


r r 

m 


(»2i': 


r      .        r  .        r      .        r  r 

cosa^siny      +sina      sin    3      cosy 
m  'm  m  m  'm 
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(a-,)*^  =  sin  a'^  sin  y^   -  cos  a'^  sin  B^   cos  y 
31m       m      m        m      m 


<»12'1 


r  r 

•cos    6      sin   Y 
tn  'm 


t^22': 


cos3      COSY      -sin   a      sinB      siny 
m  m  ni  ra  'm 


'^32'; 


r      r        r   .    r   .    r 
sin  a   cos  Y   +  COS  a   sin  6   sin  y 
m     '  m        m      m     '  m 


(^13': 


sin  3 


m 


(aoo)„  =  -sin  a  cos  6 
23'm        m      m 


(a--)!"  =  cos  a^   cos  B„ 
33  m       m      m 


The  transformation  m''  is  orthogonal  such  that 

m        ^ 


.m 


r     m      m 


(11-26) 


The  inverse  transformation  gives  us  a  means  to  express 
the  inverse  specification  set  {a.B.ylJI  with  the  original 
reference  link  r  moving  and  link  m  fixed  as  a  new  reference 
Equating  terms  from  Eq.  (11-26), 
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m 


(^la^r   =    (^Sl^m 


'23 


m        / 


'33      r 


'32 
'33  /   m 


(11-27) 


m 


'12 


^21 


n     r 


y'Ul 


m 


which    reduce    to 


3^  =  sin       (sin  a„  sin  y„  -  cos  a„  sin  B     cos  y    ) 
r  m  '  m  m  m  m ' 


m   .  -1 
a^  =  tan 
r 


m 


sin  a„  cos  y„  -  cos  a  sin  3''  cos  y'' 
-cos  oL   cos  b""        r 


(11-23) 


yJ  =  tan-^ 


cos  a„  sin  y„  +  sin  a'   sin  g""  cos  y"^ 

0] m  m m jn 

-cos  bI  cos  yI 
m    'm 


It  is  also  valuable  to  have  the  set  {a,B,Y}"  when 
shell  n  is  fixed.   The  set  can  be  obtained  by  equating  com- 
ponents of  the  matrices 
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m    r  m 


(11-29) 


to  obtain 


C  =  =-"'  t('n'r  (^sC  *  (^12)"  (=23'™  *  (^I3'"r  (^aa'l  ^ 


a„  =  tan 
m 


'^21 '"  '^13'1  *   '^22'"  '^23'!  ^  '^23'"  '"33'm 
(^3,)"r  (^13):  *  ('32)?  (^23':  ^  (^33'?  ('33': 


(11-30) 


y"  =  tan""" 
m 


(^l'"r  (=12):  *  (^12'r  (=22'^  ^3'"  ('32' 


hl'"rhl':^(^12)?(^2,':*(=13)?(^3l' 


This  transformation  allows  the  reference  link 
to  be  temporari ly  changed  so  that  the  appropriate  motion 
specification  set  {a,B ,y}  „  can  be  fed  to  the  curvature 
transformation. 

Curvature  Transformation 


A  detailed  description  of  the  3,  4  and  5  multiply 
separated  position  (MSP)  transformation  is  given  in  ref- 
erences [19,  20].   In  any  case,  the  question,  "What  is  the 
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locus  of  points  (x.y,?)  in  the  moving  shell  E  that  have 
three  positions  on  circles  about  center  (Q^  ,  Q-,  ,  Qq)  in 
the  reference  shell  E?"  is  common  and  becomes  the  central 
question  of  synthesis.   The  well-known  constraint  equation 
of  the  plane  becomes  on  the  sphere 

G(X,Y,Z)  =  02^^  +  Qi  ^i  "-   ^0   ^l  ^   ^  (11-31) 


where  D^  =  -G^,  x  +  G^^  y  -  G^^  z 


h   -    ^BZ""    '   ^7Z  y    -    ^5l' 


(11-32) 


^0    =  G-„  X  -  6,„  y  +  G.„  z 


Z  31 


4£  ^  6£ 


The  G    coefficients  of  Eq.  (11-26)  are  given  by 


Gq^  =  1  -  cos  a^cos  Q^ 


1,.  =  sin  a.  cos  y ^+    cos  a£Sin  Bg^sin  y^  (11-33) 


12 £  =  cos  ajj^sin  6^  cos  y^  -    sin  a^j^sin  Yjj, 


*Note  that  in  the  previous  references  Ani£'s  were 
used  as  G[y,£ '  s .   The  later  notation  was  adopted  to  dis- 
tinguish the  spherical  coefficients  from  the  planar  co- 
efficients. 


I  / 
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'Si   =  "^  ^"^  ^£-  ■• 


G^^  =  cos  B^sin  y^ 


^5SL   =  "'■"  V°^  ^£ 


Gg^  =  sin  6^ 


G^P  =  cos  a^cos  Y£  -  sin  a^sin  3^5 in  Y£  -  1 
Ggn  =  cos  a^sin  y^+  sin  a^sin  gj^^cos  y^ 

Equation  (11-33)  is  listed  for  finitely  separated 
positions  (FSP)  only  and  must  be  differentiated  with  respect 
to  the  independent  parameter  y   to  reflect  the  role  of  in- 
finitesimal ly  separated  positions  (ISP).   The  6  ,  notation 

ni  x» 

allows  the  same  body  of  analytics  (and  computational  tools) 
to  solve  Eq.  (11-31)  for  all  cases  of  MSP. 

There  are  four  distinctions  between  the  curvature 
transformation  program  listed  i  n  the  Appendix  and  that  cov- 
ered in  reference  [19].   First,  the  original  work  was  cen- 
tered upon  the  special  reference  transformation.   This  is 
an  effort  to  reduce  the  number  of  parameters  for  special 
cases  and  required  that  {a.,Q,y}      ^   =    {0,0,0}.   The  general 
reference  transformation  allows  the  arbitrary  choice  of 
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the  initial  position  parameters,  and  significantly  re- 
duces the  amount  of  work  required  by  the  program. 
Second,  the  homogeneous  cubic  cone 


^14  4    ^    (^24  "l  ^  ^34)  4    ^    (^44  ^1  ^  ^54  ^1  ^  ^64^  ^2 


■^  (^4  ^1  '    ^84  4    '   ^94  ^1  ^  ^104^-"°   (^^-3^) 


is  the  locus  of  homogeneous  circle  points  (z,,Zp)  in  terms 
of  known  coefficents  a..  [19].   If  (x,y)  is  a  known  point 
on  the  curve  such  as  Prj-i  »  then  Eq.  (11-34)  can  be  reduced 
to  a  quadratic.   Let 


z,  =  s  +  X 


Z2  =  t  +  y 


(11-35) 


and  substitute  into  Eq.  (11-34).   This  substitution  yields 
a  binomial  in  s,t  whose  constant  terms  collect  to  zero. 
Dropping  the  j=4  subscript  for  brevity. 


a^t  +  a2St  +  a^s  t  +  a^s^  +  a^3t^y  +  a2(2sty  +  xt^)  +  a^(2stx  +  ys^) 


+  a73s  X  +  a^t  +  a^st  +  aj^s^  +  aT3ty^  +  a,(sy^  +  2txy) 


.g. 


'1 


2  2 

+  a^(tx  +  2sxy)  +  a73x  s  +  a32ty  +  a5(sy  +  tx)  +  ag2sx  +  a^t 


+  ags  =  0 


(11-36) 


47 


Reduce  to  a  quadratic  in  polar  form  by  the  substitutions 
s  =  r  cos  e  and  t  =  r  sin  9.   Collect  terms  in  the  form 


Ar  +  Br  +  C  =  0 


(11-37) 


where 


•3         2  2  3 

A  =  a,sin  e  +  a^sin  6  cos  9  +  a^sin  9  cos  9  +  a^  cos  8 


B  =  sin  9(3aTy  +  a,x  +  a,)  +  sin  e  cos  Q{2s.^y   +  2a-x  +.ac.) 


+  cos  e(3ayX  +  a.y  +  ag) 


C  = 


2  2 

sin  9(3aiy  + 2a2xy  +  a^x  +  2a2y  +  a^x  +  ag) 


2  2 

+  cos  9(3ayX  +  2a«xy  +  a^y  +  ^agX  +  a^y  +  ag) 


and   0  <  9  <  TT 


Eq .  (11-37)  can  now  be  solved  for  r.  and  substituted  into 


the  homogeneous  coordinates 


z,  =  r.  cos  9.  +  x;  z^   =   r.  sin  9.  +  y   (11-38) 


The  G^^  coefficients  (Eq.  11-33)  are  functions  of 
Y  alone  since 


a  =  f(Y)  and  6  =  f(Y) 


(11-39) 
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The  higher  order  properties  of  the  coefficients  are  found 
by  differentiating  with  respect  to  y.      Now  suppose  the 
motion  of  a  system  is  to  be  coordinated  with  a  crank  para- 
meter (}).   This  coordination  involves  three  functional  con- 
straints 


a  =  f((})),  B  =  f((f)),  Y  =  f((}>) 


(11-40) 


For  ISP  it  is  necessary  that  the  motion  be  specified  in 
terms  of  higher  order  properties  of  these  functions  with 
respect  to  <j) .   A  new  class  of  problems  evolves  when  ^    is 
given  as  a  general  time  function 

*  =  fd) 
so  that 


a  ^  f(T),  B  =  fd).  Y  =  f(T) 


(11-41) 


Now  Eq.  (11-34)  can  represent  the  motion  specification 
set  in  coordination  with  a  general  input.   The  problem  is 
to  specify  the  time  states  of  a  shell,  and  then  synthesize 
a  system  to  meet  those  specifications.  This  is  not  a 
normal  dynamic  response  question,  since  there  is  no  pre- 
existing mechanism  responding  to  an  input;  rather  the  con- 
straints of  the  mechanism  are  sought  from  the  expanded 
curvature  transformation.   This  task  is  achieved  by 
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differentiating  the  G   's  (Eq.  11-33)  with  respect  to  time 
in  terms  of  a,  3,  y,    a,  3,   etc.   To  avoid  the  tedious 
task  of  differentiating  Eq  .  (1 1-33)  wi  th  respect  to  time, 
we  can  take  advantage  of  the  chain  rule  and  the  geometric 
derivations  of  Eq,  (11-33)  as  shown  below: 


jk 

n,^    =   — 1^  (e)  Y   =    f(T) 

dx 


T]q    =    e  =  f(Y) 


Hi    =    e   Y 


11*2    .       ,  •• 

Ho    =    e      Y      +    e   Y 


n3    =    e"'Y      +    3e"  yY    +    e'y 

A  O  O 

n^    =    e""Y      +    6e"'Y   Y    +    3e"Y     +4£"yy     +    e'T 

The  prime  superscript  on  e  denotes  derivatives  with  respect 
to  Y-   If  we  allow  e  to  represent  G   (y),  a(Y)>  or  6(y). 

Ill  X> 

then  n  represents  T   (t),  a(T),  or  B(t),  respecti vely .   If 
the  independent  parameter  is  y  3s  in  the  original  synthesis 
problem,  then 
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dx 

^- 

dy 

0 

(11-42) 


{k|k  >  1} 


One  computational  tool  can  then  be  used  for  either  geometry 
or  time  dependent  problems. 

The  tables  of  solvable  problems  include  a  column 
representing  MSP  (multiply  separated  position)  synthesis 
and  are  representati ve  of  the  possibility  of  higher  order 
geometric  specification.   This  is  in  reference  to  well-known 
structural  curvature  theory.   It  has  a  major  weakness  of 
not  being  able  to  coordinate  nonadjacent  planes  for  higher 
order  properties.   This  class  of  problems  is  solvable  only 
by  Time  State  Synthesis  (TSS). 

The  fourth  distinction  mentioned  in  the  section 
labeled  Synthesis  Tools  transforms  the  existing  FORTRAN 
program  to  a  very   efficient  interactive  APL  function.   This 
feature  greatly  enhances  the  designer's  ability  to  stay  in 
control  of  the  design  process  and  leads  to  an  earlier  solu- 
tion. 


CHAPTER  III 
DESIGN  PROCEDURES 

The  objective  of  these  design  procedures  is  to 
present  how  the  spherical  four-,  five-,  and  six-bar  mecha- 
nisms can  be  determined  dimensional ly  by  use  of  the  synthe- 
sis tools  of  Chapter  II.   There  are  a  few  hundred  potential 
design  procedures  based  on  the  problem  charts  in  Chapter 
II.   This  presentation  will  not  be  all  inclusive  but  will 
present  at  least  one  solution  for  each  problem  type.   A 
similar  effort  by  Prof.  D.  Tesar  for  planar  mechanisms 
has  proven  to  be  useful  in  both  undergraduate  and  graduate 
study  as  well  as  in  applied  industrial  practice.   The 
consistent  and  compact  nomenclature  allows  the  designer 
to  approach  complex  problems  with  fundamental  tools. 

Optimization  has  become  the  central  question  of 
planar  kinematics.   It  will  also  become  the  key  problem 
to  spherical  kinematics  as  well  as  all  fields  of  machine 
design.   Optimization  must  apply  weighting  factors  of  de- 
sign characteristics  to  a  large  collection  of  solutions. 
The  size  of  this  collection  can  be  realized  by  considering 
how  many  sets  of  motion  specifications  {a.B.y}  ai^e  used  in 
the  curvature  transformations  and  how  many  I   +    1  positions 
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are  specified.   If  three  positions  of  the  moving  plane  Z 
are  specified,  the  circle  point  locus  (locus  of  moving 
pin  joints)  is  the  entire  moving  shell.   This  shell  can 
be  uniquely  represented  by  a  pair  of  homogeneous  coordi- 
nates.  Each  coordinate  has  a  single  infinity  of  choice. 
Each  point  in  the  homogeneous  plane  maps  one-to-one  to  a 
center  point  in  the  reference  shell  E.   Each  pair  of  circle 
points  and  center  points  construct  rigid  cranks.   For  the 
three-position  problem  there  are  a  total  of  °°      cranks. 
The  locus  of  homogeneous  circle  points  for  four  positions 
of  Z  is  a  cubic  curve  (Eq.  11-34).   There  is  a  single 
infinity  of  points  on  this  curve  and  hence  there  is  a 
total  of  °°  cranks.   The  five-position  circle  point  locus 
is  found  by  the  intersection  of  two  cubics.   There  is  a 
total  of  6  cranks  for  the  five-position  problem. 

The  designer  is  expected  to  have  started  by  match- 
ing his  problem  with  the  problem  charts  of  Chapter  II. 
After  selecting  an  appropriate  mechanism  based  on  the 
number  of  free  optimization  parameters    p.,  the  nomen- 
clature  of  the  chosen  chain  should  be  referenced.   The 
procedures  then  outline  a  step  by  step  method  of  solution. 

All  of  the  procedures  will  be  produced  for  zero 
order  (displacement)  positions.   Higher  order  properties 
are  obtainable  for  each  case  by  differentiating  with  re- 
spect to  the  generic  independent  variable  x,  i.e.. 


53 


{a,e,Y} 


nk 


dx' 


{a,B,Y} 


m 


(III-l) 


^k-Jl 


where  t  can  represent  y,<I>.  or  t.   The  required  differenti- 
ation is  taken  automatically  in  the  accompanying  interactive 
computer  program.   This  statement  is  valid  for  geometric 
states  or  time  states.   It  will  be  seen  that  all  problems 
are  applicable  to  the  time  state,  but  others  may  not  be 
specified  in  the  higher  order  geometric  state.   The  objec- 
tive of  these  design  procedures  is  that  they  be  accessible 
to  the  designer  without  any  further  information  than  that 
given  in  this  document. 

1 .   Path  Synthesis 

For  path  synthesis  there  is  only  one  (F=l)  func- 
tional constraint  provided  by  the  given  path.   The  motion 
of  a  point  on  a  path  is  independent  of  a  rotation  parameter 
y.      The  standard  curvature  transformation  requires  the  set 
{a,&,y}.      The  problem  statement  gives  us 


"^  =  "e 


h  -   ^ 


Y.  =  arbitrary 


(111-2) 


This  problem  can  be  treated  by  the  colaminar  formulation 
listed  in  procedure  B  by  taking  y  as  an  arbitrary  set  of 
values . 
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2.   Colaminar  Synthesis 

This  problem  is  the  most  fundamental  and  the  most 
widely  used.   Many  of  the  complex  design  procedures  will 
reference  this  problem  but  in  no  way  detract  from  the  sig- 
nificance of  its  own  value. 

Some  point  e  of  the  moving  system  E(u,v)  moves  on 
a  path  f(U,V)=0  in  positions  e,  with  position  coordinates 
(Ug^.Vg^)  (Fig.  Ill-la).   The  transformation  between  the 
rectangular  coordinates  (X,Y,Z)  and  the  curvilinear  coordi 
nates  (U,V)  is  given  by 


sin  V 


Y  =  -sin  U  cos  V 


(III-3)* 


Z  =  cos  U  cos  V 

The  motion  of  point  e  along  the  path  is  simultan- 
eously coordinated  with  the  angular  rotation  parameter 

^^Iz'      ^"^^  ^^^  ^^eH'    ^e£'  *2£^  ^' ^  ^  ^^^  °^  Euler  angles  con- 
sistent with  the  motion  specification  set.   Hence  the  motion 
specification  set**is obtained  directly  by  the  relations: 


*Note  that  in  the  plane  Eq.  (III-3)  reduces  to 
X=U,  Y-V,  Z=l .   This  agrees  with  the  previous  planar  con- 
cepts. Equation  (III-3)  is  solved  for  the  reference  system, 
but  is  also  valid  for  the  moving  system. 

**A11  of  the  illustrations  in  this  chapter  use  the 
line  codes:   solid  link  line-specified  link;  long  dash- 
specified  shell  motion;  short  dash-constraint  found  from 
curvature  transformation. 
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(a) 


(b) 


Figure    III-l 
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"ii  =  hr  h  -   ^er  Y^ 


'2Z 


(III-4) 


It  is  now  possible  to  outline  the  procedures  for  the  basic 
mechanisms. 

2. 1   The  Four-Bar 

This  problem  depends  directly  upon  the  curvature 
transformation  to  determine  the  dimensional  constraints 
AB  and  CD  (Fig.  Ill-lb).  . 

1.  Apply  the  motion  specification  set  {U  „, 
V  .,  4^25,^   ^°  ^^^   curvature  trans- 
formation. 

2.  Select  pin  joints  B  and  C  from  the  tabulated 
ci  rcle  points. 

3.  Alternately  to  step  2,  select  pivots  A  and  D 
from  the  center  point  tabulation. 

4.  Selection  of  either  a  center  point  or  a  circle 
point  determines  the  other.   Points  B  and  C 
are  both  in  the  rigid  moving  shell  and  can 

be  joined  by  a  rigid  l"'nk  of  arc  length  BC. 
It  is  identical  to  make  the  rigid  link  360°- 
BC.   The  pivots  A  and  D  are  both  points  of 
the  rigid  reference  E  and  provide  the  base 
link  of  the  four-bar  ABCD. 

3.   Angular  Coordination  of  Two-Link  Chain 

The  system  of  two  links,  2,3  moving  relative  to 
link  T  (Fig.  III-2a)  is  completely  known  if  the  following 
data  are  given 


^'  '^h'  ^3£ 


£  =  0,1 ,2, 


(III-5) 
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u 


(b) 


Figure    III-2 
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The  curvilinear  coordinates  for  the  state  of  b  in  z  are 
given  by  Eqs.  (II-8  and  9),  where  e-.  =  ^pa' 


a 


^'-    i^ 


tan 


1  /  sinabsinii^p 
cos  aB 


2A   H.k 

UT 


sin"   (sin  a  b  sin  ^7.    ) 


The  rotation  of  link  3  about  pivot  b  completes  the  motion 
specification  of  link  '2. 


dx 


^3  *  P2 


where  q^q    i ^  given  by  Eq.  (11-12)  as 


'21         .    k 
dx 


tan"   (tan  y]^2     c°s  ^^) 


3.1   The  Four-Bar 


The  four-bar  can  be  used  if  the  symbolic 
transfer  (Fig.  III-2b) 

a->A,  b^B,  T-^1  ,  2-^2,  3^3 


occurs 


59 


2.   The  curvature  transformation  can  be  used 
to  determine  crank  CD. 


4.   Angular  Coordination  of  Two  Cranks 

The  problem  of  "function  generation"  deals  with 
synchronizing  the  angular  motion  of  two  cranks  rotating 
about  fixed  pivots  a  and  b  (Fig.  III-3a).   All  of  the  mecha- 
nisms under  study  are  capable  of  this  performance.   We 
shall  limit  the  procedures  to  the  four-bar  and  five-bar 
1.   It  is  felt  that  these  have  enough  versatility  to  ac- 
complish most  angular  motion  problems  of  this  type. 

4.1   The  Four-Bar  by  Inversion 


In  Fig.  III-2  we  have  the  motion  of  2  and  3 
given  relative  to  T.   If  we  jxow  invert  the 
£urrent  statement  by  fixing  2  and  allowing 
1  to  move  we  make  the  symbolic  transfer  (Fig 
III-3b). 


a-»-A 
1^1 


b^d 

Ul  ,  3^4 


The  moti^on  specification  for  the  two-link 
chain  (1  and  3)  is  given  by 


{a, 8,7)4^  ^  Eqs. (11-8,9,10,12) 

3.  From  the  curvature  transformation  we  find 
circle  point  C  and  center  point  B  to  form 
the  rigid  constraint  of  link  3. 

4.  Invert  the  mechanism  back  to  the  original 
system  by  fixing  link  1. 

Note  that  the  initial  positions  of  ^1.   and  il/r  are  not  met 

by  arcs  AB  and  CD.   However,  the  finite  displacement  values 

Ai|)2   and   L^-^   are   met   by   moving    shells    2   and    4. 
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(a) 


U 


a  =  7r-v//^ 


(b) 


Figure    1 1 1-3 
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4. 2  The  Five-Bar  1  by  Inversion 


The  required  symbolic  transfer  is 

a->A,    b->F,    c->-D 

T^l ,    2^2 ,3^5 

asinFig.III-4. 

The  overall  gear  train  value  from  A  through 
F  to  D  is  Q=MN.   Links  2  and  4  are  coupled 
by  this  train.   The  influence  of  e-.  and  e„ 
on   e^  can  be  seen  by 

(e 


3Z 


^30^  =  Q(9U-  ^lo)  "^(Q2£-Q20)     (^"-6) 


We  have  no  control  on  030  so  that  we  may  ar- 
bitrarily set  it  to  zero.   Eq.  (II-6)then 
yields  83^. 

The  motion  specification  for  link  4  becomes 

{a,&,y}l^   ^  Eqs.  (11-14,16) 

The  curvature  transformation  gives  pin  joint 
C  as  the  circle  point  locus  and  point  B  as 
tjie  center  point  locus.   The  pair  form  arc  ' 
BC  and  provide  the  needed  constraint  for  the 
f  i  V  e  -  b  a  r  1  . 

Invert  the  mechanism  by  fixing  link  1. 


Path  Coordination 


The  four-bar  has  enough  parameters  to  treat  the 
coordination  of  a  crank  with  a  path  (Fig.  III-5a).   Planar 
four-bars  use  the  path  cognate  tool  that  is  not  available 
on  the  sphere.   Without  this  tool  we  suffer  the  weakness 
of  using  six-bars  to  solve  a  problem  that  a  four-bar 
should  be  able  to  treat.   The  procedure  for  the  six-bars 
follows: 
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Figure    III-4 
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U 


(a) 


U 


(b) 


Figure    III-5 
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1.   Choose  ab  and  find  the  position  of  b  by 

hz  =   ^^'^  ^u  '^'^  ^^ 


Y.  „  =  -cos  Q,^    sin  ab 


(III-7) 


I.„    =   cos  ab 
bit 


From  curvilinear  coordinates  of  e(Ug,Vg) 
find  the  position  e(X  ,Y  ,Z  )  by 


\l   -   -^i"  ^;i   ^°^    ^eii 


3, 
4, 


•eJi 


"^   ^£   "^    ^eil 


From   dyad    analysis    (Chapter    IV),    find    c(X,Y,Z), 


^1%   ^"^  ^Sil 


(Fig.  IIl-5b),  using  dyad  bee 


Use  any  of  the  procedures  for  the  coordination 
of  a  three-link  chain  to  follow. 


6.   Angular  Coordination  of  a  Three-Link  Chain 


The  coordination  of  a  three-link  chain  is  the  first 
problem  of  this  presentation  that  cannot  treat  higher  or- 
der geometric  properties  by  direct  use  of  the  S   's  in  the 
curvature  transform.   The  coordination  of  two  nonadjacent 
shells  for  higher  order  position  properties  can  only  be 
treated  by  considering  another  independent  parameter  other 
than  e^^  =  Y^  (Fig.  II-6). 

It  is  more  useful  in  this  problem  to  consider  the 
crank  angle  ti^pj  =  Qig  =  4)  as  independent  and  specify  all 
other  necessary  position  functions  in  terms  of  this  parameter 
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(b) 


Figure    II 1-5 
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If  ({)  =  f(T),  the  motion  is  time  dependent  and  the  spheri- 
cal T   's  must  be  used  in  terms  of  the  specification: 
mil  ^ 

The  state  of  shell  3    can  be  expressed  from  Eqs.  (11-8,9,10) 


a 


hi         .  k 
di 


tan 


sin  a  b  sin  0 ,  \ 


cos  ab 


/ 


'b£    ,  k 

di 


sin"    ( s  i  n  a  b  s  i  n  0 ,  ) 


b£    ,  k 
dx 


'2   ^  Pb 


^^^•"^   Pb£  =  TT 
dx 


tan'  (tan  6,   cos  ab) 


The  state  of  shell  4  can  be  expressed  from  Eqs.  (11-14,16), 


a 


cl         ,  k 
dx 


tan-^  [  -h    '^"^  ^1 


Y^  cos  6, 


cJi    .  k 
dx 


sin"    (X^  cos  e^   -  Y2„  sin  6,  ) 


c£    .  k 
dx 


'3   +  °c   ^  Pc 
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where        a 


c£         ,    k 

UT 


tan 


,/x,^ 


clx 


tan"    (sin    3        cos    a      ) 


6. 1   The  Stephenson  3,  Case  1 


Coordinate  the  problem  statement  and  the 
Stephenson  3  by  the  symbolic  transfers 
(Fig.  III-6b) 


a^-F,    b-^-D,    c->C 
T-^3,    2-^1,    3^4, 


4^5 


Find  the  first  motion  specification  set  for 
link  4  as 


3,  From  the_curvature  transform  find  the  con- 
straint GD  (Fig.  III-6b). 

4.  Find  the  second  motion  specification  set  for 
link  5  as 


{a.B.Ylg^  =  ^^c£'^r^cil^ 


5.  From  the_curvature  transform  find  the  con- 
straint BA  (Fig.  III-6b)  . 

6.  Pin  joints  C,D,  and  G  form  a  rigid  ternary 
Pin  joints  B  and  C  form  a  rigid  binary  to 
complete  the  Stephenson  3  mechanism. 
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6 . 2   The  Stephenson  3,  Case  2 


Make  the  symbolic  transfer  (Fig.  III-7a) 

a-^-A,    b-»-B,    c^-C 

U3,    2^2,    3^5,    4^4 

Find  the  motion  specification  set  for  link 
4  as 


From  the  curvature  transform  find  two  con- 
straints GH  and  DF.   Pin  joints  C,  G,  and  D 
form  a  rigid  ternary  supported  by  four-bar 
HGDF  and  driven  by  binary  AB  (Fig.  Iir-7a). 


6  .  3  The  Stephenson  2 


1.  Make  the  symbolic  transfer  (Fig.  III-7b) 

a-^A,    b-^F,    c^D 

T^2,    2-v3,    3^1,    4-v4 

2.  Find  the  motion  specification  set 

{a,6,Y}4^  =  ^°'c£'^cr^c£^ 

3.  From  the_curvatu re  transform  find  the  con- 
straint CB. 

4.  Find  the  motion  specification  set  for  link  3 


(a.B.Ylgj^  =  (0,0, e^^} 


Find  the  inverse  motion  specification  set 
when  link  3  is  fixed  as 


{a,3,Y}2j^  =  ^°'°'-^U^ 
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(a) 


U 


(b) 


Fi  gure    1 1 1-7 


>0 


Find  the  motion  specification  set  for  link  4 
wi  th  link  3  f i  xed  by 


{a.e.y} 


4£ 


fn  =  3 
Eq.  (11-30)  where  {r   =    2 

!m  =  4 


From  the  curvature  transform  with  the  specifi- 
cation set  from  step  6,  find  link  6  from  circle 
points  G  and  center  points  H. 

Invert  the  mechanism  to  the  original  system 
by  fixing  link  2.   Note  that  points  A,  F, 
and  H  form  a  rigid  ternary  as  well  as  points 
C ,  D ,  and  G. 


6 . 4  The  Stephenson  1 


Make  the  following  symbolic  transfer  (Fig 
III-8a) 


a->F, 
T^l, 


b-^A, 
2^3, 


c^B 
3^2, 


4^5 


Find  the  motion  specification  set  of  link  5 
with  reference  link  1. 


From  the  curvature  transform  find  link  4  with 
circle  point  C  and  center  point  D. 

The  coordinates  B(X,Y,Z)  can  be  found  from  sub- 
stituting a^^  and  Bg.  into 


hi   =    ^■'^  ^5£ 


'bsL    =  ^■'■"  «5£  ^°^  hi 


hz    =  "^  ^5ii  ^°^  ^5£ 
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(a) 


(b) 


Figure    III-8 
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The  coordinates  of  C(X,Y,Z)  in  the  £  +  1  po- 
sitions can  be  determined  from  dyad  analysis 
using  dyad  BCD.   Analysis  will  also  provide 
the  relative  angle  ^  FDC  (6^). 

The  problem  now  leaves  the  determination  of 
link  6  to  the  procedure  of  function  genera- 
tion between  links  3  and  4  given  in  Section 
4. 


6.5  The  Watt  1 


1.  Make  the  symbolic  transfers  (Fig.  III-8b) 

a^A,    B^D,    c->H 

T^l  ,    2-^2,    3^3,    4->5 

2.  Find  the  motion  specification  set  of  link  4 
by 

3.  From  the  curvature  transform  find  constraint 
link  3  with  circle  point  C  and  center  point 
B. 

4.  The  motion  specification  set  of  link  2  is 
given  by 

{a,B,Y}  2^  =    {0,0,0-,^} 

5.  The  motion  specification  of  link  5  is  given 
by 

{a.B,Y}^£  =  ^"c£'^c£'^c£^ 

6.  Invert  the  mechanism  by  fixing  link  2. 

{a,3,Y}^£  =  {0,0, -e^^} 
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The  motion  of  link  5  with  respect  to  link  2 
is  given  by 

Cn    =    2 


{a.B.Y}^^  ->  Eq.  (11-30)  where 


r  =  1 
m  =  5 


8.  From  the  curvature  transform  find  the  con- 
straint link  6  with  circle  point  F  and  cen- 
ter point  G. 

9.  Points  A,  D,  and  G  construct  a  rigid  ternary 
link  (2).   Points  C,  D  and  H  construct  a  rigid 
ternary  link  (4).   Inverting  to  the  original 
system  by  fixing  link  1  satisfies  the  problem. 


7 .   Angular  Coordination  of  a  Crank  with  a  Two-Link  Chain 

The  motion  of  a  two-link  chain  whose  angular  state 
is  coordinated  with  a  crank  is  closely  associated  with  the 
three-link  chain.   The  problem  is  depicted  in  Fig.  III-9a. 
The  basic  specification  follows  that  of  the  three-link 
chain  (cab  -t->  abc  ) . 

7 . 1   The  Stephenson  2 


1.   Make  the  symbolic  transfer  (Fig.  III-9b) 


a^B,    b->C,    c^A 
U2  ,    T-^5  ,     3-^4  , 


4^3 


Fix  link  3  by  allowing  Q-\i    =    tt  -^21   ^^^  2£  " 
^31   to  find  the  motion  specification  for  link 
4. 

{a, 3, 7)4^  ->  {a^^.B^^.Y^^} 

From  the  curvature  transform  find  constraint 
links  DF  and  GH.   Points  D  and  G  come  from 
the  circle  point  laws  and  points  F  and  H  come 
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(a) 


^z-n 


(b) 


Figure    II 1-9 
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from  the  center  point  locus.   Points  C,  D 
and  G  construct  the  rigid  ternary  link  4. 

Invert  the  mechanism  to  the  original  system 
by  fixing  link  2. 


7.2  The  Stephenson  1 


1.   Make  the  symbolic  transfer  (Fig.  Ill-lOa) 


a-»-D, 
T-1, 


b-»-C,    c^F 
2-^4,    3^5, 


4-^3 


2.   Allow  the  rotation  parameters  of  the  three- 
linkchaintobe 


®u  =  ^  -  ni 


2£ 


=    -n    +   i\) 


2£ 


'3£    ^3£ 


and  consider  its  inversion  by  fixing  link  3. 

3.   Find  the  motion  specification  set  for  link  4 
relative  to  link  3  as 


{a,6,Y}4^  =  ^°'a£'^a£'Ya£^ 


4.  From  the  curvature  transform  find  link  6  with 
circle  point  G  and  center  point  H.   Note  that 
points  C,  D,  and  G  form  the  rigid  ternary  link 
4. 

5.  Find  the  motion  specification  set  for  link  5, 


{a,6,Y}5^  =  ^"b£'^b£'^b£^ 
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Figure  III-IO 
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From  the  curvature  transform  find  link  2  with 
circle  points  B  and  center  points  A.   Note 
that  the  points  H,  A,  and  F  form  the  rigid 
ternary  1 i  nk  3. 

Invert  the  mechanism  back  to  the  original  sys- 
tem by  fixing  link  1 .  . 


8.   The  Angular  Coordination  of  Three  Cranks 

-  This  problem  can  be  considered  as  the  double  func- 
tion generator  and  can  be  solved  directly  from  applying 
the  four-bar  function  generator  procedures  twice.   Only 
three-pivot  mechanisms  are  eligible  for  this  problem. 

8.1   The  Watt  2 

1.   Make  the  symbolic  transfers  (Fig.  Ill-lOb) 


a^A,    b^-D,    c->G 
T-^2  ,    2-^1  ,    3->4 , 


4->6 


Consider  the  problem  parameters  ij/o^  and  if^^g, 

and  use  the  coordination  of  two  cranks  in 

Section  4;  find  link  Swith  pin  joints  C  and 
B. 

Consider  the  problem  parameters  ip^i   and  ij^^j^ 
and  use  the  coordination  of  two  cranks  in 
Section  4  to  find  link  5  with  pin  joints  F 
and  H  . 

Points  D,  C,  and  H  are  all  on  rigid  ternary 
link  4  to  couple  the  two  four-bars  together. 


The  Angular  Coordination  of  Three  Adjacent  Planes 


If  link  4  is  considered  as  a  reference,  the  prob- 
lem is  transferred  to  the  type  of  the  three-link  chain  with 
(bac  -!->  abc)  . 
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9.1   The  Watt  1 


1.   Make  the  symbolic  transfers  (Fig.  Ill-lla) 

a->A,    b^-D,    c^G 

T->1  ,    2^2,    3^6,    4  =  4 


Let    e-|ji   =    ^2i   ^"^    02£   "^    ?4Jl   t°    find    the   motion 
specification    of    link    4. 


1 


{a,3,Y}4ji    =    ^^biL'^b£'^b£^ 


From  the  curvature  transform  find  link  3  with 
circle  point  C  and  center  point  B. 


4.   Let 


'U 


TT  -  ? 


4ii 


02,  =  TT  +  T 


'3£  ~  ^32, 


and  find  the  motion  specification  set  of  link 
6  with  respect  to  link  4  (considering  link  4 
as  fixed) 


{a,B,Y} 


en 


^""cZ'^Ci'^Cl^ 


From  the  curvature  transform  find  link  5  with 
circle  point  F  and  center  point  H. 

Note  that  points  C,  D  and  H  form  the  rigid 
ternary  1  ink  4  . 


9.2   The  Watt  2 


1.   Make  the  symbolic  transfer  (Fig.  Ill-llb) 
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U 


(b) 


Figure    1 1 1- 11 


ao 


3, 
4 


a^D,    b^H,    c-^-C 

T^2,    2->4,    3-^5,    4  =  3 


Let  e-jn  =  4)2£'  ^21      ^    ?4il  ^^^    '•'''"'^  ^^^  motion 
specifications  for  link  3. 

From  the  curvature  transform  find  the  constraint 
link  1  with  circle  point  B  and  center  point  A. 

Let  Q-[i   =   i|j^£  +  1^  and  625,  =  .C5£  and  find  the 
motion  specification  set  for  link  5. 

{a, 3,7)5^  =  {"cr^cr^cji^ 


5.   From  the  curvature  transform  find  the  constraint 
link  6  with  circle  point  F  and  center  point  G. 

10.   Double  Colaminar  Synthesis  with  Common  Knee  Joint 

This  problem  has  had  several  valuable  applications 
in  planar  kinematics.  It  is  expected  that  it  will  also  be 
valuable  on  the  sphere. 

10.1   The  Stephenson  3 

1.  Make  the  symbolic  transfer  (Fig.  III-12a) 

e->C,    3->4,    2->5,    T-^3 

2.  Defind  the  motion  specification  set  forlink 
5  as 


{a,8,Y}^,  =  ^U^'^c'^P^ 


^1 


5  ^B 
/2 


I  / 


(a) 


(b) 


Figure    III-12 
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3.  From  the  curvature  transform  find  the  constraint 
link  2  with  circle  points  B  and  center  points  A. 

4.  Define  the  motion  specification  set  for  link  4 
as 


ia,B,y} 


4Jl 


=  ^'ch^nh 


From  the  curvature  transform  find  the  two  con- 
straint links  1  and  6.   Note  that  points  D, 
G,  and  C  form  ternary  link  4.   The  resultant 
mechanism  is  depicted  in  Fig.  III-12a. 

Switching  the  roles  of  wy   and  ^J   by  the  sym- 
bolic transfer  2->4  and  3->5,  the  steps  2  through 
6  can  be  repeated  to  double  the  number  of  solu- 
tions (Fig.  III-12b). 


11 .   Colaminar  Motion  Coordinated  with  a  Crank 

This  problem  involves  the  coordination  of  a  moving 
shell  with  the  motion  of  a  crank.   It  may  be  the  most  prac- 
tically useful  procedure  of  this  treatise. 

11.1   The  Stephenson  2 

1.  Make  the  symbolic  transfer  (Fig.  III-13) 

a->A,     e^D 

T^2 ,     2->3 ,     3-^4 

2.  Define  the  motion  specification  set  for  link 
4  as 


{ct,3,Y}4j^  "  ^^D'^D''^3^ 


From  the  curvature  transform  find  the  con- 
straint link  5  with  circle  point  C  and  center 
p  0  i  n  t  B  . 


d3 


U 


Figure    1 1 1- 13 
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From  analysis  of  dyad  DCB  find  the  motion  of 
link  5  and  angular  motion  <i)c^.  Also  find  the 
base  angle  4)2  and  angles  4)4, 

Let 


U 


ni^  *2£ 


^21  ^  ha 


'3Z 


'u 


and  find  the  motion  of  link  4  with  respect  to 
link  3  by  considering  link  3  fixed. 


{a,6.Y}4^  =    ^"c£'^c£'^cil^ 


6.  From  the  curvature_transform  fi_nd  the  system 
constraints  links  DP  as  1  and  GH  as  6. 

7.  Invert  the  mechanism  back  to  the  original  ref- 
erence of  link  2. 


CHAPTER  IV 
SPHERICAL  DYAD  ANALYSIS 

Kinematics  is  based  upon  developing  an  under- 
standing of  the  motion, independent  of  time, of  intercon- 
nected rigid  bodies.   An  analysis  technique  is  necessary 
to  describe  the  motion  of  any  point,  shell,  or  relative 
angle  in  a  mechanism.   This  chapter  will  generate  a  collec- 
tion of  tools  to  provide  a  concise  and  efficient  means  of 
interpreting  the  analysis  problem.   It  should  be  noted  that 
analysis  is  a  necessary  compliment  to  synthesis  and  design, 
and  may  well  be  the  final  critical  step  of  optimization. 

The  use  of  dyads  for  analysis  of  planar  mechanisms 
has  proven  to  be  the  most  efficient  and  versatile  closed 
form  tool  that  is  presently  available.   Referring  to 
Fig.  IV-1,  a  dyad  is  represented  as  a  pair  of  rigid  links 
(BC.CD)  connected  by  a  common  joint  at  point  C.   If  the 
time  states  of  the  two  remaining  ends  (B,D)  are  known, 
then  the  time  state  of  C  can  be  found.   Note  that  BCD 
forms  a  spherical  triangle.   Without  loss  of  generality, 
we  will  consider  only  unit  spheres  such  that  6,c,  and  8  are 
all  unit  position  vectors  of  points  B,C,  and  D  respectively, 
The  arc  BD  can  be  found  by  the  scalar  product  of  6  and  3. 
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Figure    I V -  1 
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BD  =  cos"\b-d) 


(IV-l) 


The  interior  angle  ^  CBD  (0)  can  be  written  in  terms  of 
the  spherical  cosine  lav; 

cos  CD  =  cos  BC  cos  BD  +  sin  BC  sin  BD  cos  0.     (IV-2) 

Solving  for  cos  6,  and  taking  into  account  the  possible 
range,  we  find 


-1  < 


cos  6 


COS  CD  -  cos  BC  cos  BD 
sin  BC  sin  BD 


<  1 


(IV-3) 


Satisfaction  of  this  condition  defines  whether  or  not  the 
dyad  can  be  physically  connected. 

In  order  to  find  the  position  vector  c,  we  must 
first  rotate  6  to  (0,0,1)  by  a  rotation  u  about  the 
-X-axis  and  then  a  rotation  v  about  the  -y-axis.   This 
transformation  can  be  described  in  matrix  form  by 


0 
0 
1 


M^B 


(IV-4) 


where 
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COS  V 

sin  u  sin  V 

M,  = 

0 

cos  u 

sin  V 

-sin  u  cos  V 

-cos  u  sin  V 

sin  u 
cos  u  cos  V 


However,    from   the    curvilinear    to    rectangular    transforma- 
tion , 


X    =    si n    V 
y   =    -sin    u    cos    v 
z    =   cos    u    cos    V 
Rd  =  Vl-x2    =    CO 


(IV-5) 


^B 


S  V 


Substituting  Eq.  (IV-5)  into  Eq.  (IV-4),  the  new  position 
of  D  is  found  by 


R 


B 


"^B^B 


■^B^B 


(IV-6) 


Now  rotate  d '  to  lie  in  the  {-y  ^z)    plane  by  a  rotation  tf) 
about  the  -z-axis. 


0 

Zq., 


M^a* 


(IV-7) 


where 
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M, 


cos  4" 

sin  (() 

0 

-sin  (j) 

cos  (t> 

0 

0 

0 

1 

However  the  position  of  8'  can  be  defined  by 


x^r  =  sin  BD  sin  ^ 
ypj,  =-sin  BD  cos  ^ 
Zq,  =  cos  BD 


(IV^8) 


and 


sin    BD    =yi-Zn.       =    VxlT+yT^ 


D 


'D 


D 


Rp, 


Substituting  Eq.  (IV-8)  into  Eq.  (IV-7), 


M„  = 


2   R 


-yQi 

'^D- 

0 

1 

Rq. 

-Xp, 

-^D' 

0 

0 

0 

Rq. 

(IV-9) 


The  position  of  c"  is  given  by  the  elementary 
expressions , 


Xp„  =  sin  BC  sin 


y-,,  =-sin  BC  cos 


C" 


(IV-10) 


Zp„  =  cos  BC 
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Both  transformation  matrices  M,  and  M,,  are  orthogonal 

I      c 

such  that  their  inverse  is  identical  to  their  transpose 
{m7  =  M-|  and  mZ  =  Mp)   Reversing  through  the  inverse 
rotations,  we  find  the  position  vector  c  of  point  c  in 
the  original  system  position  by 


c 


m]  M^  c" 


(iv^n) 


The  assumptions  of  rigidity  of  links  BC  and  CD, 
as  well  as  the  constraint  of  the  unit  sphere,  impose  the 
following  set  of  conditions  on  points,  B,C,  and  D. 


b*c 


c  'd  = 

b.b  = 

As  /\ 

c»c  = 


d-d  = 


cos  BC  =  constant 

cos  CD  =  constant 

cos  0  =  1 

cos  0=1 

cos  0=1 


(IV-12) 


Differentiating  Eq.  (IV-12)  yields  a  set  of  con- 

•   •  • 

ditions  for  the  velocity  vectors  6,c,  and  3. 


G*c  +  G-c  =  0 

c'3  +  c-3  =  0 

btb  =  0 


C'C  =  0 


(iV-13) 


3.3  =  0 
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These  can  be  collected  in  several  ways  to  avoid 
singularity  for  the  solution  of  c,  e.g., 


^B    ^B    ^B 


^D    -^D    "D 


J 


^B    ^B    ^B 


>^D    ^D    ^D 


M  *-  -^ 


(IV-14) 


If  the  determinant  of  M(det  M)  is  nonzero  then  premulti- 
plying  Eq.  (IV-14)  by  -M   is  sufficient  to  find 
(Xp.ypjZp).    If  det  M  is  zero,  then  a  different  collec- 
tion of  Eq.  (IV-13)  should  be  assembled  in  Eq.  (IV-14). 
To  find  the  acceleration,  we  can  differentiate 
Eq.  (IV-13)  again  to  give  the  set. 


b'C  +  2  b«c  +  b'C 


/\        /v. 


d*c  +  2  d*c  +  d«c  = 
c  •  c  +  c  •  c  = 


0 
0 
0 


(IV-15) 


b«b  +  b«b  =  0 

•    ■       •< 

d'd  +  d-d  =  0 


These  two  can  be  selected  in  several  combinations 
to  f ind  c ,  e.g . , 


i)2 


^B    ^B 

■^b" 

'^c] 

-2xg 

'h 

^h~\ 

r-    t 

^D     ^D 

^"d 

^c 

+ 

2xo 

^h 

2^D 

^c 

0      0 

0 

.^c. 

_    \ 

k 

^C- 

-^C 

^B  ^B  ^B 


>^D  ^D  ^D 


\     ^C  ^CJ 


^C 


L^CJ 


M 

(IV-16) 

If  M  in  Eq.  (IV-14)  is  n on  singular,  then  Eq.  (IV-16)  can 
be  premul  ti  pi  ied  by  -M"   to  find  (xp,yr.,Zp).   This  pro- 
cedure could  be  carried  further  to  find  higher  order 
properties  if  desired. 

Now  that  we  have  shown  how  to  find  the  time  state 
of  point  C  in  the  dyad  with  points  B  and  D  given,  an  ob- 
servation can  be  noted  on  their  extension  to  points  all  on 
the  same  rigid  shell.   If  this  is  the  case,  such  that  BD 
and  e  are  constant,  the  same  procedure  can  be  utilized  to 
find  the  time  state  for  a  third  point  C  on  the  same  body. 
This  means  that  one  tool  can  treat  both  dyads  of  two  rigid 
links  or  ternary  links  of  one  rigid  shell. 


Prismatic  Equivalence 

The  slider  or  prismatic  joint  is  relatively  common 
in  planar  mechanisms.  A  slider  is  generally  considered  to 
be  a  link  that  is  moving  along  a  linear  path  with  a  center 
of  rotation  at  infinity.   Links  in  the  plane  transform  into 
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great  circles  on  the  sphere.   By  rotating  a  90°  crank  on 
the  sphere,  any  great  circle  of  that  sphere  can  be  gen- 
erated.  Hence,  any  90°  crank  on  the  sphere  is  the 
equivalent  of  a  slider  in  the  plane. 

The  general  spherical  dyad  enables  us  to  illus- 
trate and  discuss  the  four  four-bar  chains  (RRRR,  RRRP, 
RRPP,  and  RPRP)  and  all  of  their  inversions  by  the  use  of 
the  spherical  equivalents  to  the  five  planar  dyads 
RRR,  RRP,  RPR,  PRP,  PPR). 

Four-Link  Chain  Equivalence 


The  RRRR  chain  is  the  general  four-bar  linkage. 
Both  the  planar  and  spherical  representations  are  given 
in  Fig.  IV-2a.   All  inversions  are  the  same  and  all 
spherical  links  are  different  from  90°. 

The  RRRP  chain  is  shown  in  Fig.  IV-2b.   One 
link  (CD)  is  inf ini te  in  the  planar  case.   The  offset  in 
the  planar  chain  from  pivot  A  to  straight  line  £,  is 
h  =  AA' .   The  corresponding  offset  in  the  spherical  chain 
is  A  =  AA'.   In  the  planar  case  AM)  =  CD  =  «>.   In  the 
spherical  case  A'D  =  CD  =  90°.   No  slider  needs  to  be 
represented  physically  at  C  although  C  always  travels  on 
a  great  circle  c-j  .   The  travel  on  Z^    is  measured  by  S^  . 
The  corresponding  measure  in  the  spherical  case  is  a. 
Note  the  link  AD  =  h  +  0°  in  the  planar  case,  and 
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Figure    IV-2 
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AD  =  A  +  90°  in  the  spherical  case.   There  are  two  unique 
inversions  of  the  RRRP  chain.   The  first  occurs  when  link  1 
(or  link  4)  is  fixed.   It  is  called  the  slider-crank 
mechanism  (Fig,  IV-3a).   The  second  inversion  (Fig.  IV-3b) 
occurs  when  you  fix  link  2  (or  link  3). 

The  third  common  four-link  chain  is  the  RRPP 
(Fig.  IV-4a),  which  has  three  unique  inversions.   Note 
that  link  1  contains  straight  line  constraints  z-i    and  £_ 
in  the  planar  case.   This  means  that  in  the  spherical 
configuration  z^    and  z^    become  great  circles  c,  and  Co 
which  have  a  constant  relative  orientation  angle  e  . 
The  corresponding  angle  in  the  plane  is  0p.   Great  cir- 
cles have  their  centers  at  the  center  of  the  sphere,  and 
normals  to  the  planes  of  the  great  circles  intersect  the 
sphere  at  points  A  for  c-j  and  point  D  for  c^.   Since  these 
points  (D,A)  correspond  to  c-j  and  c^,  both  in  link  1, 
points  D  and  A  are  in  link  1.   Hence  DA  is  a  constant  arc 
length  equal  to  the  orientation  angle  e  .   Fixing  link  1 
makes  the  double-slider  mechanism  which  generates  ellipti- 
cal motion  on  the  plane.   Fixing  link  2  (or  link  4)  forms 
the  Scotch  Yoke  Mechanism  which  is  well  known  as  a  harmonic 
input  output  function  generation.   Fixing  link  3  makes  the 
double  oscillating  block  mechanism  and  is  the  last  of  the 
unique  inversions  for  this  chain.   In  this  case,  pin 
joints  B  and  C  become  fixed  pivots.   All  points  in  link  1 
generate  cardiod  motion. 
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Figure    IV- 3 
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A  at  CD 


(a)     RRPP     Four    Link    Chain 


A  at  CO 
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Figure    IV-4 
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Figure    IV-5 
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The  last  four-link  chain  is  the  RPRP  chain  (Fig 
IV-6).   It  has  only  one  unique  inversion  which  is  the 
slide-guide  mechanism.   In  this  case,  there  are  two 
straight  line  constraints  £-,  and  %^    each  of  which  may 
be  considered  in  each  of  the  four  links  of  the  chain. 


Dyad  Collection 

The  RRR  dyad  (Fig.  IV-7a)  is  the  generic  dyad  dis- 
cussed in  the  earlier  sections.   In  both  the  plane  and  the 
sphere  it  is  constructed  of  rigid  links  1  and  2  with 
revolute  joints  B,C,  and  D. 

The  RPR  dyad  (Fig.  IV-7b)  has  joint  C  at  infinity 
in  the  plane.   There  is  one  linear  constraint,  straight 
line  i^i  ,  in  link  1  with  offset  h-j  from  B.   This  corre- 
sponds with  great  circle  c,  in  link  1  for  the  spherical 
case.   The  offset  A,  is  measured  from  pin  joint  B.   Note 
that  point  C  could  be  taken  on  either  side  of  the  great 
ci  rcl e . 

The  RRP  dyad  (Fig.  IV-7c)  has  joint  D  at  infinity 
in  the  planar  configuration.   There  is  one  linear  con- 
straint l^    offset  from  pin  joint  F  by  h2  =  FF'.   In  the 
spherical  configuration,  the  line  l^    ^"^    represented  by 
great  circle  c^    in  link  2.   Point  F  is  offset  by  Ap  =  FF ' 
from  C2  and  CD  =  90°. 
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Figure    IV- 7 
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The  PRP  dyad  (Fig.  IV-8a)  has  both  pin  joints  B 
and  D  at  infinity  in  the  planar  case.   This  configuration 
has  similar  characteristics  to  the  RRP  dyad. 

The  RPR  dyad  has  two  linear  constraints.   In  the 
plane,  straight  lines  £■,  and  £„  can  be  considered  to  lie 
in  the  same  link,  in  this  case  link  1.   Also  pin  joints 
B  and  C  lie  at  infinity,  perpendicular  to  lines  £,  and 
^2*    respectively.   Because  of  this,  pin  joints  B  and  C 
may  be  considered  to  both  lie  in  link  1  (as  always  has 
been  the  case  for  all  dyads).   Point  A  is  offset  by  dis- 
tance h  =  AA'  from  line  £, .   The  lines  £,  ,£p  correspond 
to  great  circles  c,  ,Cp  in  link  1  of  the  spherical  repre- 
sentation.  Offset  AA '  becomes  A-,,  and  the  dyad  becomes 
BCD  where  AB  =  A^  +90°,  and  CD  =  90°. 

As  is  apparent  in  all  cases  of  the  spherical 
dyads,  BCD,  all  joints  are  pin  joints.   The  only  distinc- 
tion occurs  when  90°  is  the  arc  length  between  these  pin 
joints.   Hence  only  one  analysis  routine  is  sufficient 
to  interpret  the  very  large  class  of  planar  and  spherical 
mechanisms  that  can  be  constructed  by  dyads.   To  illustrate 
an  example  of  an  assemblage  of  dyads,  consider  the  eight- 
link  mechanism  of  Fig.  IV-9.   The  planar  configuration  is 
built  from  four  dyads, 

1.  BCD  (RRR) 

2.  BEC  (RPR) 

3.  EFG  (PPR) 

4.  HJK  (RRP) 
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Figure    IV-8 
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H(R) 
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EC=90*+ A2 


Figure    IV- 9 
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The  spherical  configuration  is  built  from  four  RRR  dyads 
with  the  following  link  lengths. 


1. 
2. 
3, 
4, 


BCD   BC,  CD 

BEC   90°  +  A, ,  90°  +  A, 


'1 


EFG   e^  ,  90° 
HJK   HJ,  90° 


If  the  arc  lengths  AB  ,  BC  ,0) ,  A-,  ,  A2  .HJ  are  quite 
small  (largest  ~.5°)  then  a  planar  approximation  with 
errors  less  than  10'  can  be  expected.   The  error  can  be 
further  refined  by  reduced  magnitudes  on  the  largest  link 
length  that  is  not  90°. 


Spherical  Six- Bars  by  Dyads 

Figure  IV-10  illustrates  a  collection  of  common 
mechanisms  that  can  be  solved  directly  by  dyads.  Each 
mechanism  is  built  from  a  four-bar  substructure  and  then 
combinations  of  dyads  applied  to  the  ternary  points  e- 
and  to  fixed  pivot  g  to  provide  the  analysis  of  a  large 
percentage  of  mechanisms.   The  primary  exception  to  this 
analysis  procedure  is  the  Stephenson  2  mechanism. 

If  oj^  (or  ojg)  is  the  input  to  the  Stephenson  2 
shown  in  Fig.  IV-11,  then  the  input  loop  is  based  upon  a 
pentagon  rather  than  the  quadrilateral  loop  common  to  all 
other  six-bar  mechanisms.   This  means  that  the  Stephenson  2 
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Figure    IV- 10 
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cannot  be  built  from  dyads  directly,  but  is  dependent  upon 
the  higher  level  Assur  group  shown  in  Fig.  IV-12.   No 
closed  form  analysis  tool  is  available  for  this  class, 
and  in  order  to  avoid  iterative  techniques,  an  alterna- 
tive closed  form  solution  will  be  developed. 

Consider  the  Stephenson  3  inversion  by  fixing  link 
3  and  allowing  link  2  to  be  free  to  rotate  about  point  A 
as  shown  in  Fig.  IV-11.   If  we  now  drive  link  1,  the 
Stephenson  3  mechanism  can  be  completely  analyzed  by  the 
use  of  dyads.   Note  in  the  Stephenson  2,  that  point  B  is 
a  fixed  pivot  and  has  no  motion.   If  a  transformation 
matrix  could  operate  on  point  B  of  the  Stephenson  3  to 
rotate  it  back  to  Bg,  then  this  same  transformation  could 
be  applied  to  points,  C,D,F,6,  and  H  to  give  their  inverted 
motion  in  the  Stephenson  2.   This  tactic  is  valid  for  posi- 
tion and  all  higher  order  properties  as  well.   Assume, 
without  loss  of  generality  that  point  A  is  at  (0,0,1),  then 
the  rotation  matrix  from  B  to  Bq  is  given  by 


6q  =  M  6 


(IV-17) 


where 


M  = 


"'  \  *'>'  \\\'^-\  "^ 


\'^ 


Xk  Yu  '  X,   X.  +y.    y^ 
bg  •'b  I  b   bg  ■'b   •'t 

1 


tR 
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Figure    IV-11 
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Figure  IV-12 


no 


where   R      =vl-zf     =   constant 

Differentiating    Eq.    (11-17) 


bQ=0=Mb+Mb 


b=o=Mb+2Mb+Mb 


ClV-18) 


This  technique  is  supported  by  an  APL  function  listed  in 
the  Appendix,  and  has  all  of  the  planar  implications  that 
apply  to  previous  tools. 

Body  Rotations 

If  the  time  states  of  two  points  on  a  rigid  body 
are  known,  the  velocity  and  acceleration  poles  can  be 
determined.   The  velocity  pole  is  the  instantaneous  axis 
about  which  pure  rotation  exists.   Referring  to  Fig. 
IV-13,  this  rotation  axis  ui  can  be  related  to  the  transla 
tion  of  a  point  on  the  body  by  the  vector  cross  product. 


a  =  u)  X  a 


(IV-19) 


This  can  be  broken  into  three  components  and  then  expressed 
as  a  matrix  of  simultaneous  equations 


y 

2  J 


-X 


A 


u. 
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•-  z  -' 


(IV-20) 
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Matrix  A  is  singular  and  cannot  be  used  alone  to 
solve  for  uj.   Writing  Eq.  (IV-20)  for  each  of  the  two 
given  points  A,B  gives  six  equations  from  which  three 
could  be  chosen  to  give  a  nonsingular  matrix  A.   Two  more 
equations  can  be  added  by  noting  that  the  linear  and  angu- 
lar velocities  are  perpendicular  such  that 


a  •  0)  =  0 


(IV-21) 


For  computational  purposes,  the  singular  condition  can 
be  most  easily  avoided  by  summing  three  sets  of  three 
equations ,  e.g.. 
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(IV-22) 


^C  ^C  \ 
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L  (jJ ) 


The  angular  acceleration  u  can  be  readily  found  by  differ- 
entiating Eq.  (IV-22)  with  respect  to  time. 

Planar  motion  possesses  the  characteristic  of  equal 

angular  velocities  and  accelerations  for  all  lines  through 

any  point  in  a  rigid  body.   In  Fig.  IV-9,  if  w  and  Z   are 


113 


found  for  the  pole  axis,  then  y   and  y  about  point  A  have 
the  same  direction  and  magnitude  as  w  and  u  respectively. 
This  is  due  to  the  fact  that  the  direction  cosine  between 
parallel  axes  is  unity.   However,  the  sphere  does  not 
share  this  characteri stic,  and  the  absolute  rotation  proper- 
ties about  a  point  other  than  the  pole  must  account  for 
the  direction  cosines  of  the  angle  between  the  position 
vectors  of  point  A  and  the  pole  axis  (Fig.  IV-14). 

Relative  Angles 

In  planar  and  spherical  mechanisms,  the  relative 
angle  between  adjacent  links  is  often  of  significant 
interest.  This  interest  may  be  due  to  the  consideration 
of  transmission  angles,  or  for  function  generation.  Con- 
sider in  Fig.  IV-14,  the  relative  displacement  angle  be- 
tween links  BC  and  CD  about  point  C.  Writing  the  cosine 
law  for  the  spherical  triangle  BCD, 


cos  BD  =  cos  BC  cos  CD  -  sin  BC  sin  CD  cos  9-    (IV-23) 


Solving  for  0q  we  have. 


ep=  cos 


cos  BC  cos  CD  -  cos  BD 
sin  BC  sin  CD 


(IV-24) 
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y  =  Q*(jj  =  leal  cos  A| 
y  =  d'CJ  -\(jj\  cos  X2 


Figure    I V-1 4 
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where  BC  and  CD  are  given  and  BD  can  be  found  by  Eq. 
(IV-1).   The  result  of  Eq.  (IV-24)  is  an  angle  between  0 
and  IT,  but  in  reality  9^  may  range  from  -tt   to  ir.   The 
appropriate  sign  can  be  found  from  the  sign  of  c  •  (8x8), 
Note  that  only  the  position  of  point  C  is  needed  to  find 
the  proper  sign  of  6q,  and  no  higher  order  properties 
are  used.   Differentiating  Eq.  (IV-24)  with  respect  to 
time  yields. 


BD  sin  BD 


sin  dr.   sin  BC  sin  CD 


(IV-25) 


Differentiating  Eq.  (IV-1)  and  substituting  into  Eq.  (IV- 
25)  yields 


^0  = 


6'3  +  b'3 


sin  6^  sin  BC  sin  CD 


(IV-26) 


Differentiating  Eq.  (IV-25)  once  more  yields  the  relative 
angular  acceleration. 


b»d  +  2b'd  +  b'd 
sin  Gp  sin  BC  sin  CD 


-   0p  cot  Qr. 


(IV-27) 


This  process  could  be  continued  without  end  as  long  as  the 
higher  order  properties  of  points  B  and  D  are  known. 
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System  Inputs 

A  simple  crank  (Fig.  IV-15a)  is  the  most  common 
input  to  a  linkage  system.   Given  the  location  of  pivot 
A  and  the  arc  length  AB ,  the  time  state  of  point  B  can  be 
found  for  the  chosen  values  of  6  while  being  driven  at  a 
constant  velocity  m.      The  coordinates  of  B'  for  A'  = 
(0,0,1 )  are  given  by 


Xni    =  sin  e  sin  AB 

D 


yn. 


■COS  e  sin  AB 


(IV-28) 


Zn.  =  COS  AB 


The  rotation  matrix  M-,  in  Eq.  (IV-4)  can  be  used  to  rotate 
B'  to  the  original  system  position  B.  Let  the  B  subscript 
of  Eq,  (IV-6)  be  replaced  by  A  and  find. 


6  =  MJ  B' 


(IV-29) 


Equation  (IV-29)  can  be  differentiated  to  yield  the  higher 
order  properties  of  point  B.   If  AB  is  equal  to  90°,  then 
the  same  formulation  describes  the  slider  input  on  a 
great  circle. 

The  cycloidal  crank  (Fig.  IV-15b)  offers  a  second 
type  of  input  to  the  dyad  system.   Given  the  location  of 
A  with  links  AA '  and  A'B  and  gear  ratio  N,  the  coordinates 
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Figure    IV- 15 
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of   coordinates    of   B'    with    A   temporarily    held    at    (0,0,1) 
are    given    by 


Xg,   =  X2  cos  e  -  Yp  sin 


yg,   =  X2  sin  e  +  Y2  cos 
2g.    =  \ 


(IV-30) 


where 


X„  =  sin  A'B  sin  A 


Y2  =  -(sin  AA'   cos  A'B  +  cos  AA'   sin  A'B  cos  X) 
Z2  =  cos  AA'   cos  A^  -  sin  AA'   sin  A^  cos  A 


(IV-31) 


A  =  Aq  +  N  (6-6^). 

The  pivot  may  be  moved  to  a  general  location  A,  by 
Eq.  (IV-29). 

All  of  the  above  tools  have  been  collected  into  a 
group  of  compact  and  efficient  APL  functions  to  provide 
interactive  access  to  the  spherical  and  planar  analysis 
problems.   This  set  of  tools  create  the  desired  output 
for  all  input  positions  simultaneously.   The  use  and  list- 
ing of  these  programs  are  described  in  the  Appendix. 


CHAPTER  V 
DISPLACEMENT  ANALYSIS 

The  previous  chapter  discussed  methods  of  analyz- 
ing a  specified  closure  branches  of  a  mechanism  including 
those  with  inputs  that  cannot  rotate  360°  due  to  closure 
limitations.   A  higher  level  of  understanding  can  be 
achieved  by  treating  algebraically  all  angular  input- 
output  relations  of  a  given  mechanism  simultaneously. 

This  section  will  seek  the  development  of  input- 
output  functions  as  well  as  all  significant  relative 
angles  of  six-link  spherical  mechanisms.   A  similar  treat- 
ment for  four-  and  five-link  mechanisms  is  fully  treated 
in  references  [21  ,  23] . 

Watt  1 


For  the  Watt  1  mechanism  (Fig.  V-1),  the  cosine 
law  for  quadrilateral  1237  yields. 


^12  ~-    '°^  °^37 


(V-1) 


where 
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Figure    V-1 
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Z,5  =  sin  a-jT  (X,   sin  e,  +  Y,   cos  Go)  +  cos  a-,o  Z 


•12 


*23   '1   "■"  °2       '1 


23  "1 


X-.   =  sin  ay-j   sin  0, 


Y-,   =  -(sin  a-ip  cos  a^,   +  cos  a-.^  sin  a-,-.   cos9,) 


Z-j   =  cos  a-.2  cos  a^,   -  sin  a-.^  sin  a,,   cos  9, 


If   e,    is    known   as    an    input   then   X,,    Y,  ,    and    Z,    are   all 
known.      Substituting    the    half-tangent    identities 


2x. 
sin  0.  =  — ~    ;  cos  0- 
1+x^  ^ 


1^ 


(V-2) 


(where  x-  =  tan  -^) 


into    Eq.    (V-1),    terms    can    be   collected    in   a    quadratic 
input-output   equation    of   the    form 


Ax2   +    Bxp    +   C 


(V-3) 


where 


A  =  cos  a23  Z^   -sin  a^o  Y,   -cos  a^. 


B  =  2  sin  ttp.  X 


C  =  cos  a^^  Z^   +  sin  a^,  Y^   -cos  a.., 
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Let  9o   =   ■^3  +  63*    tlien    the   half-tangent   laws   yield 


^3  = 


^'^"  "37  '  ^12 
^12 


^12 


sin  a37  +  Y^^ 


(V-4) 


where 


X,2  =  X-,  cos  62  -  Y"!   sin  62 


Y,o  =  cos  a^,   (X,   sin  e^  +  Y-,   cos  So)  -sin  a^o     Z 


12 


^23    VA-,    :>...  ^2  -    .^ 


23     "1 


The   angle   87   can    be    found    by   another    set   of    half-tangent 
laws    of   the    form 


^7  = 


S"*"  ^37  ~  ^21 
^21 


'21 


sin  a37  +  Y2^ 


(V-5) 


where 


Xpi   =  X2  cos  9,   -  T2  sin  9, 


Y2^   =  cos  a^^    (X^  sin  9-,  +  Y2  cos  6^)  -sin  c^^■^   1^ 


X2     =  sin  a^o  sin  e^ 


Yo     =  -(sin  a, 2  cos  a23  +  cos  a, 2  cos  a^^  cos  92) 


Zp     =  cos  a, 2  cos  Up^  -sin  a, 2  sin  a23  cos  62 
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In  Fig.  V-2,  let  e'7  =  ^,-9-,  and  consider  the  penta- 
gon formed  by  34567  and  write  the  cosine  law 


Z734  =  cos  agg 


(V-6) 


where  1^^^   =    sin  a^^    (X73  sin  e^  +  Y73  cos  e^)  +  cos  a^^    Z^. 


X73  =  X7  cos  63  -  Y^  sin  63 


'73  "  "^  "34  (^7  s"""  63  +  Y7  cos  0")  -  sin  a   Z 


^73  "  ^^'"  ^^34  (^7  ^■'"  83  +  Y-,  cos  e^)  +  cos  a3^  Z^ 


X7   =  sin  ag^  sin 


Y7   =  -(sin  a37  cos  a^^  +  cos  a37  sin  ag^  cos  e^) 

Zy   =  cos  037  cos  ag7  -  sin  a37  sin  a^^  cos  e^ 

With  X73,  Y^3,  and  2^3  known,  Eq.  V-6  becomes  a  quadratic 
in  x^  with  coefficients. 


^  ^  ^°'  "45  ^73  -  "^'^  "45  Y73  -  cos  a 
B  =  2  sin  a^,  X 


56 


45   73 


(V-7) 


C  =  cos  a^g  Z73  +  sin  a^g  Y73  -  cos  a^g 
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Figure    V-2 
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Finally  angles  0^,  and  Gg  +  Tg  can  be  found  by  the  half- 
tangent  laws 


^5  = 


^734 


734 


^^■"  ^'56  '    ^34 


(V-8) 


where  X^^^   =    X^^  ^°s  84  "  Y73   sin  64 


734 


=  cos  a^g  (X73  sin  e^  +  Y^^  ^os  64)  -  sin  a^g  Ij. 


and 


X  ' 


^"■"^  ^56  -  %37 
^437 


437 


^^■"  ^56  '   ^37 


(V-9) 


where  X^^y  =  ^43  *=°s  e^  =  Y^^  sin  Q' 


Y-07  =  cos  a^-,  ( X.o  sin  e4  +  Y 


437 


67  '''43  -^  ■"  °7    '43 


cos  6^)  -  sin  ag7  2^3 


^43  =  ^4  <^°s  63  -  Y^  sin  93 


^43  "  ^°^  «37  ^^4  sin  63  +  Y^  cos  63)  -  sin  a37  Z 


37  -4 


2^3  =  sin  037  (X4  sin  63  +  Y^  cos  95)  +  cos  a37  2^ 


X^   =  sin  a3.  sin  9. 
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Y^  =  -(sin  a.g  cos  a^^  +  cos  a.^  sin  a^^  cos  6^) 


Z4  =  cos  a^g  cos  034  -  sin  a^^  sin  a^^    cos  e^ 


x^  =  tan  --^- 


Note  that  since  quadratics  were  solved  in  x^  and  x.,  there 
are  four  possible  closures  for  an  input  9,. 

Watt  2 

The  Watt  2  mechanism  can  be  broken  into  two 
four-bars  of  quadrilateral  loops  I  and  II  as  in  Fig.  V-3. 
Write  the  cosine  law  for  loop  I, 


^17  ~-   '°'  ^6 


(V-10) 


where  Z,^  =  sin  a-,n    (X-]  sin  9'  +  Y-j  cos  9^)  +  cos  a^^  Z. 


X"!   =  sin  agi  sin  e^ 


Y"!  =  -(sin  a-.-,    cos  ag-,  +  cos  a-jy  sin  agi  cos  Q') 


Z^  =  cos  a-^-j    cos  agi  -  sin  a^y  sin  ag^   cos  9]) 
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Figure   V-3 
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e-  =  8^  +  -1 


Write  the  cosine  law  for  loop  II, 


^27  =  "'  ^34 


(V-11) 


wh 


ere  1^-^    =  sin  a-^  (X^  sin  e!|  +  Y2  cos  9")  +  cos  a,^  Z 


SI n  apo  SI n 


Y2   =  -(sin  a27  cos  a^,  ■*"  cos  0.^-.    sin  a^o  cos  6') 


Z2   =  cos  a27  cos  a23  -  sin  0.^-^    sin  a23  cos  6' 


®2   "^  ^2  "^  "^2 


If  e-i  is  the  input,  Eq.  (V-10)  can  be  written  as 
a  quadratic  in  x^  with  coefficients 


A,  =  cos  a  g  Z-j  -  sin  a  j.  Y-.  -  cos  a^ 


B^  =  2  sin  ay n  X, 


(V-12) 


C^  =  cos  a^g  Z^  +  sin  a^^  Y^  -  cos  a^g 
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Substituting  e"  =  t^  +  t^  -  6'  into  Eq.  (V-ll) 


gives 


cos  a^^  =  sin  a,y  {X2(sin(T^+T,)cos  0'-cos(T^+T-)sin  9') 


+  Y2(cos(Ty+Ty)cos  0'+sin(T^+T^)sin  6^)}+  cos  a.y  Z_ 


(V-13) 


Equation    (V-3)    can    be   written   as    a    quadratic    in 


Xy   with   coefficients, 


Ag  =  cos  ct  y  1^  -  cos  a^^  -  sin  ct-y(X     sin(Ty+T-,)+  "i^  cosd^+Ty)) 


B2  =  2  sin  a. 7(72  sinCxy+x,)-  t^   cosCxy+Ty)) 


(V-14) 


C2  =  cos  a^y  Z2  -  cos  a^^  +  sin  a-^CX  sinCx^+x^)*  Y2  cosCx^+Xy)) 


Eliminating  x'  between  quadratic  Eqs.  (V-12)  and 
(V-14)  yields  a  quadratic  in  Xp  by  the  expansion  of  the 
resultant, 


^    ^1   ^    ^1 


Mq       0       0       o 

^    ^1    ^1    ^1 


=  0 


(V-15) 


M(^      ^O        O        O 
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Stephenson  1 

The  Stephenson  1  input-output  equation  could  be  found 
in  at  least  two  different  ways.   The  input-output  is 
based  upon  the  quadrilateral  1276  and  could  readily  be 
written  by  introducing  the  constant  ternary  angles  at 
vertices  1  and  2.   However,  it  is  felt  that  including  as 
many  sides  as  possible  gives  a  more  pallatable  result. 
It  should  be  recognized  that  the  numeric  results  must  be 
identical . 

Consider  the  hexagon  123765  of  Fig.  V-4  and  write 
the  cosine  law. 


^5123  =  "'  %?  ^^-^^^ 


where  l^^^^    =    sin  a^^i)^^^^    '"'"  ^3  '  ^12  '°'  ^3^^  ^°'°'73  ^512 


^512   =  ^1  "'  ^2  -  ^1  '■■■"  ^2 


^12   =  ^°^  "23^^51  '■*"  ^2  ^    h^    "'  ^2^-    '■•■"  ^23  ^51 


^512   =  '^"  ^23(^1  ^■'■"  ^2  ""  ^1  ^°'  ^2^^  ^°'  °'23  ^51 


Substituting  Xg-io'  ^512'  ^"^  ^512  ^"^°  ^^-  (V-16) 
and  collecting  terms. 
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Figure    V-4 
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(X^i  Y3  +  Y^^  J^)   sin  02  -(Xg^  X3  -  Yg^  Y3)  cos  6^ 


•(^51  ^3  -  "'  ^67^  =  ° 


(V-17) 


where  Xo  =  sin  a^^  sin  a-. 


Y3  =  -(sin  apn  cos  a^^  -  cos  a23  sin  a^^  cos  x^) 


Z2  =  cos  a23  cos  a-j2   ■*"  sin  ap3   sin  a.^^   cos  t„ 


^sr  ^5  '°'  ^1  -  ^1  '"■■"  ^1 


Yg^=  cos  a^2(^5  ^^"  ^1  "^  ^5  ^°^  ^i^"  ^^"  "i?  ^5 


Zg^=  sin  ai2(Xc  sin  0^  +  Y   cos  e,)+  cos  a, 2  Z^ 


X5  =  sin  a^g  sin  x^ 


Yr  =  -(sin  agi  cos  a^c    -    cos  ac-,  sin  aj-,   cos  Xj.) 


Zg  =  cos  ar-j  cos  a^g  +  sin  a^i  sin  a^^    cos  Xc 


If  9-1  is  taken  to  be  a  known  input,  then  e^  is 
the  only  unknown  in  Eq.  (V-17).   Using  the  half-tangent 


133 


identities,  Eq.  (V-17)  can  be  written  as  a  quadratic 
with  coefficients. 


'^    ~-    ^51^3  "  ^1^3   -  ^51^3   ^  "'  «67 


B  =  2(X5^Y3  +  Y5^X3) 


(V-18) 


^  =  -^51^3  ^  ^51^3  ■  ^51^3  ""  ^°'  "'e? 


For  completeness,  the  other  angles  can  be  found 
by  the  appropriate  use  of  sine,  sine-cosine,  and  cosine 
laws.   Writing  the  cosine  law 


^4  =  hi 


(V-19) 


where  Z,,  =  sin  a,-(Xi  sin  9-  +  Y,  cos  e^)  +  cos  a^^  Z 


12 


t23VA-,  ....  ^2  ■  -1 


23   1 


X"!   =  sin  an •]  sin  9, 


Y^   =  -(sin  ai2  cos  a^^    +  cos  a^2  ^in  a^i   cos  9,) 


Z^   -  cos  a^2  cos  a^^  -  sin  ai2   sin  ag,   cos  9-] 


^4  "  ^°^  "^34  '^°^  °'45  ■  ^^'"  ct34  s^'"  ^45  ^os  9^ 
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After  solving  the  input-output  Eq.  (V-17),  then  Eq. 
(V-19)  can  be  written  as  a  quadratic  in  x^  with  coeffi 
c  i  e  n  t  s  , 


A  =  cos  a^^  cos  a^g  +  sin  a^^    sin  a^^  -  Z^^ 


or      A  =  cosCa^^-a^^)-  I^^ 


(V-20) 


B  =  0 


C  =  cos(x34+a45)-  Z^^ 


Writing  the  half-tangent  laws  for  9-,  Sr.  6^,  and  6-, 

o    b    D        / 


(^2^^) 


Xj2"  ^4 


(X^2  ^^4)    ^2  -  U 


(V-21) 


(Y^i  +  Y4)    X^T  -  X^ 


^^5""  (X27+  X4)  -  Y 


21  '    ^4 


(V-22) 


^123  ~  ^'''"  "'S? 
^5123 


'5123 


^123  "■  "■■■"  ^67 


(V-23) 
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^215    ~    ^'''"   ^7 
^3215 


'3215 


^3215    '    ^^"   \7 


(V-24) 


Stephenson  2 

In  the  Stephenson  2  in  Fig.  V-5  write  the  cosine 
law  for  the  pentagon  loop  12345, 


^345  =  ""^^   "12 


(V-25) 


where  Z^^^  =  sin  ^51 (X^^  sin  6^  +  Y^^  cos  9^)+  cos  a^^  Z 


X34  =  X^  cos  84  -  Y3  sin  9^ 


Y34  =  c°s  °'45^^3  ^^'"  ^4  "^  ^3  ^°^  ^4^"  ^^'"  ^45  ^3 


l-iA      =    sin  a/,(-(X,  sin  9;,  +  Y^  cos  9.  +  cos  a.c    Z 


■34 


'45^"3  "'"  ^4    '3 


45  "3 


X^   =  sin  a23  sin  9^ 


Y3   =  (sin  034  cos  a23  +  cos  a^^    sin  a^^   cos  83) 


X3   -  cos  034  cos  a23  -  sin  a3.  sin  a^^    cos  83 


136 


Figure    V-5 
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Applying  the  half-tangent  identities  to  Eq. 
(V-25),  a  quadratic  in  x^  can  be  written  with  coeffi- 
cients. 


^1  =  '°'  ''5^    ^34  -  '■■■"  Si  ^4  -  '°^  °'12 


^1  =  2  sin  a^^    X3^ 


(V-26) 


^1  =    "'  ^'Bl  ^34  ^  ^^"^  Si  ^34  -  "^  «12 


At  this  point  either  63  or  e^  may  be  chosen  as  an  input, 

but  this  would  leave  the  other,  as  well  as  6^,  unknown 

b 

and  hence  cannot  be  solved  yet. 

Writing  the  cosine  law  for  the  pentagon  loop  67345, 


^345  =  "'  ^67 


(V-27) 


where  Z^^^  =  sin  a^^U'^^    sin  <p^   +  Y'^  cos  ^^)+   cos  a^g  Z^ 


34 


X34  =  X^  cos  9^  -  Y^  sin  e^ 


Y' 


34   "  ^°^  °'45^^3  ^^'"  ^4  "^  ^3  ^°^    ^4^"  ^^'"  ^t^g  Z^ 


^34   "  ^■'"  "45^^3  ^^""  S4  +  Y'  cos  64)+  cos  a^^  Z^ 
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;*      =    sin  a^-    sin   (j)_ 


Y^      =    -{sin   a^.    cos   a^^    "^    ^°^   "34    ^''"   cty,    cos   ({>-) 


Z^      =    cos   a^^   cos   a^^    ~    ^"'"   ^1,4    sin   a^o    cos    (f)- 


and 


^    =    ^    ^    ^5 


*3    =    ^3    "   ^3 


(V-28) 


Substituting    Eq.    (V-28)    into    Eq.    (V-27)    and    apply- 
ing   the    half-tangent    identities,    a    quadratic    in   Xc    can    be 
written   with   coefficients. 


A2  =  -  cos  ag7  +  sin  a^^CX^^  sin  t^  -  Y'^  cos  t^)+  cos  a^g  Z' 


34 


B2  =  2  sin  a^gCX^^  cos  x^  -  Y^^  sin  t^) 


C2  =  -  cos  ag7  -  sin  a^^U'^^   sin  T5  -  Y^^  cos  t^)^  cos  a^g  Z^^ 


(V-29) 


139 


Equations    (V-26)    and    (V-29)    represent    a    pair   of    quadratics 
inX(-oftheform 


as 


N  X5  +  B^  Xg  +  C^  =  0 


B2  Xg  +  B2  Xg  +  C2  =  0 


(V-30) 


Eliminating  x^.  the  resultant  can  be  expressed 


A,  B, 


h    ^2 


A,  C, 


A,  C, 


0  0 


^1  ^1 


B2  C2 


=  0 


(V-31) 


Expanding  Eq.  (V-31)  yields  an  eighth  degree 
polynomial  in  x^  and  x^.   Selecting  either  as  input,  the 
other  can  then  be  solved  for, 

Stephenson  3 

For  the  Stephenson  3  in  Fig.  V-6,  write  the  cosine 
law  for  hexagon  123765, 


^3215  =  "'  ^'e? 


(V-32) 
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Figure    V-6 
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wh 


ere  732^5  =  sin  a^^{X^^^    sin  t^  -  ^ ^^^    cos  T3)+cos  a^g  Z 


321 


^321   =  ^32  "'  ^1  -  ^2  "'■"  ^1 


^321   =  "'  °'51^^32  '■■■"  ^1  ""  ^2  ^°'  ^1^-  '■•■"  "51  ^32 


•321   "  ^^"  "51^^32  ^^"  ^1  "^  ^32  ^°^  ^1^"^  ^°^  "51  ^32 


^32   ~  ^3  ''"^  ®2  ~  ^3  ^^"  ®2 


32 


cos  a,p(X  sin  6,  +  Y_  cos  6,)-  sin  a,^    I 


12 


2    "3 


12   3 


Z32   =  sin  a^2(^3  S''"  ^2  "^  ^3  ^°^  ^2^"^  ^°^  "l2  ^' 

Applying  the  half-tangent  identities  to  Eq.  (V-32),  a 
quadratic  in  x-]  can  be  written  with  coefficients 


^1  =  ^32^1  ■  ^2^  -  ^32^1  ^  ^°^  "67 


^1  =  2(^32^1  '    hzh^ 


(V-33) 


■1  ^  -^32^1  ^  ^2^  -  ^32^1  ""  ^°^  "67 


If  82  is  a  known  input,  then  0,  is  the  only 
sii 
Eq.  (V-17) 


unknown  since  X^,  Y^,  and  Z^  are  constants  included  in 
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For  the  pentagon  12345,  we  have  the  cosine  law. 


^321  =  "'  "45 


(V-34) 


where  I^^^    is  given  within  Eq.  (V-32).   Eq.  (V-34)  may 
also  be  expressed  as  a  quadratic  in  x-j  with  coefficients, 


^2    =   cos  a5^  Z32  -  sin  a^^  Y32  -  cos  a^^ 


B2  =  2  sin  ag^  X32 


(V-35) 


C2  =  cos  a^^  Z32  +  sin  a^^  Y32  -  cos  a^^ 


Eliminating  x,  from  the  quadratic  Eqs.  (V-33) 
and  (V-35)  yields  the  resultant, 


^  ^1  ^  ^^ 


rtrt   *^Q        O        ^O 


=  0 


(V-36) 


Expanding  Eq.  (V-36)  gives  a  multinomial  which  is 
eighth  degree  in  62  and  fourth  degree  in  e..   For  any 
specified  value  of  Q^,    we  obtain  a  maximum  of  four  values 
of  83  and  vice  versa.   Corresponding  values  of  9-,  can  be 
found  by 
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*1 

^1 

h 

h 

*2 

'^2 

^2 

'^2 

^ 

h 

^ 

<:i 

«2 

h 

«2 

^2 

{V-37) 


Closure 


The  degree  of  the  input-output  equation  is  a 
direct  means  of  indicating  the  maximum  number  of  real 
closures  for  the  mechanism.   If  a  mechanism  is  degree 
eight  in  the  input-output,  then  there  is  a  maximum  of 
eight  geometric  configurations  possible  in  which  to  assem- 
ble the  links.   For  the  Stephenson  1,  there  are  two  pos- 
sible reflections  (2-6,  2-5)  for  a  given  input  9 -j  .   This 
yields  a  total  of  four  configurations  or  closures  and  is 
identical  if  02  is  chosen  as  input.   The  Stephenson  2 
input-output  relation  is  based  upon  the  pentagon  loop 
34567  and  has  a  maximum  of  eight  real  configurations. 
This  was  clearly  established  earlier  in  reference  [29].  The 
Stephenson  3  closures  can  be  explained  in  a  similar  manner 
considering  the  input  loops  1276  and  73456.   If  9_  is 
considered  as  a  known  input  then  the  six-link  mechanism 
is  being  driven  by  a  quadrilateral  loop  and  has  a  maximum 
of  four  real  closures.   If  9„  is  considered  as  a  known 
input  then  the  mechanism  is  being  driven  by  a  pentagon 
loop  and  has  a  maximum  of  eight  real  closures. 
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All  of  the  above  formulation  is  readily  reducible 
to  planar  results  and  all  closure  conditions  are  directly 
anal ogous. 


Comparison  of  Chapters  IV  and  V 
Analysis  Techniques 


The  analysis  tools  of  the  preceding  chapters  have 
presented  the  designer  with  two  basic  alternatives  to 
iterative  analysis.   If  iteration  is  required  at  the  level 
of  examination  of  each  mechanism,  and  then  again  to  loop 
through  a  collection  of  mechanisms,  the  cost  and  accumula- 
tive error  can  become  impressively  large.   The  designer 
must  be  able  to  treat  his  problem  in  the  most  direct  and 
efficient  means  available.   Closed  form  analysis  is 
necessary  whenever  possible, and  iterative  solutions  should 
be  pursued  as  a  last  resort.   This  leads  to  expanded 
efforts  in  the  development  stages  but  one  should  take 
the  attitude  that  'once  it  is  done,  it's  done.' 

All  of  the  included  techniques  are  developed  with 
this  spirit.   It  may  be  helpful  to  readers  of  this  work 
to  outline  a  comparison  of  the  approaches  in  Chapters  IV 
and  V. 
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Dyad  Analysis 

1.  Based  upon  a  combination 
of  vectors  and  spherical 
trigonometry 

2.  Treats  one  closure  at  a 
time 

3.  Treats  all  states  with 
respect  to  input  posi- 
tions simultaneously 

4.  Is  easily  accessible  to 
higher  order  states 

5.  Finds  translational  and 
angular  states  in  either 
an  absolute  or  relative 
reference 

6.  No  iteration  necessary 
per  mechanism 


7.  Dead  centers  are  handled 
but  not  predicted 


Displacement  Analysis 

1.  Based  upon  spherical 
tri  gonometry 


2.  Treats  all  closures  simul- 
taneously 

3.  Treats  one  input  position 
at  a  time 


4.  Is  awkward  for  higher 
order  states 

5.  Finds  relative  angular 
displacement 


6.  Solution  of  polynomials 
(above  quartic)  is  itera- 
tive, but  polynomial 
development  is  closed  form 

7.  Easily  develops  an  under- 
standing of  special  proper- 
ties such  as  dead  center 
predictions 


The  user  must  decide  what  type  of  information  is 
necessary  in  his  examination.   Most  often  the  dyads  ap- 
proach will  prove  more  appropr  i  ate,  but  the  displacement 
analysis  can  give  a  higher  level  of  understanding  of  the 
nature  of  closure  loops. 


CHAPTER  VI 
APPLICATIONS 

Three-Link  Chain 

The  three-link  chain  is  a  valuable  aid  to  several 
synthesis  procedures.   A  complete  sample  solution  will 
follow  to  validate  the  use  and  operation  of  the  chain 
problem.   This  example  will  follow  the  design  procedures 
of  Section  6.1  in  Chapter  III.   The  understanding  of  this 
problem  has  become  the  key  to  understanding  basic  principles 
of  spherical  synthesis. 

Consider  the  three-link  chain  in  which  the  grounded 
link  is  driven  by  a  motor  at  a  constant  angular  velocity.* 
Assume  that  the  relative  angles  62  and  63  (Fig.  Vl-la) 
are  functions  of  the  input  crank.   The  desired  discrete 
values  to  be  satisfied  are  given  by 


*Inputs  are  not  restricted  to  only  constant  veloc- 
ities and  may  be  general  input  functions.   This  fact  is 
due  to  the  time  state  formulation  and  was  a  previous  re- 
striction of  geometric  synthesis. 
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Figure   V I  -  1 
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i 

k 

^U 

^2Z 

^3£ 

uni  ts 

0 

0 

30 

-5 

10 

deg 

1 

1 

120 

-50 

50 

deg/sec 

2 

2 

0 

-100 

0 

deg/sec^ 

3 

0 

60 

-25 

28 

deg 

Let  arc  aB  =  30°  and  arc  be  =  40°.   The  problem  is  now 

fully  specified.   The  state  of  shell  3  is  given  by  the 

motion  specification  set  {o..„y&.„,y,„}    as  listed  in  Sec- 

b£  bx,  bic 

tion  6  of  the  design  procedures.   This  formulation  has 
been  coded  in  the  APL  function  LINK3. 


Z  k  '^^  5M  ZM  units 

0  0  25.56  14.48  21.57               deg 

1  1  -27.71  53.67  60.85  deg/sec 

2  2  -113.94  -51.91  -46.38  deg/sec^ 

3  0  16.10  25.66  31.31      deg 


This  motion  set  can  be  treated  by  the  curvature  transfor- 
mation for  4  MSP.   This  is  automatic  in  the  LINK3  function. 
The  resultant  is  a  set  of  cranks  satisfying  the  conditions 
of  the  motion.   Recall  that  these  motion  specifications 
are  generated  from  the  O-j^  and  92£  data,  and  are  independent 
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of  635.   Selection  of  crank  GH  from  the  locus  provides 
a  four-bar  FDGH  that  is  sufficient  to  control  the  motion 
of  Q21   ^^'^^  system  input  6^^^.   From  the  locus,  circle 
point  G  and  center  point  H  are  chosen  (Fig.  Vl-lb). 

G(X.Y,Z)  =  (-.0364,  .2409,  .9699)^^q 
H(X,Y,Z)  =  (.3658,  .1211,  .9228) 


D(X,Y,Z)  =  (.25,  -.433,  .866)^^ 


F(X,Y.Z)  =  (0,0,1) 


GH  =  24.38° 


DG  =  43.41 


CG  =  83.38° 


The  state  of  shell  4  (Fig.  Vl-la)  is  given  by  the  motion 
specification  set  ^<^cl'^c!i''^cl'^    ^^    listed  in  Section  6. 
This  computation  is  also  an  integrated  part  of  LINK3. 


a 


cl 


'cl 


Y 


CS. 


units 


0 

0 

SI  .11 

24.85 

17.09 

deg 

1 

1 

-16.66 

79.19 

59.89 

deg/sec 

2 

2 

-73.55 

-76.30 

-21 .87 

deg/sec^ 

3 

0 

61  .28 

39.26 

33.97 

deg 
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Applying  this  data  set  to  the  curvature  transformation 
yields  the  locus  of  cranks  BA.  From  the  locus,  circle 
point  B  and  center  point  A  are  chosen  (Fig.  VI-2). 


B(X,Y,Z)  =  (.0564,  .2826,  .9576) 


1  =  0 


A(X.Y,Z)  =  (.0878,  -.1993,  .9760) 
C(X,Y.Z)  =  (.4203,  -.8400,  .3432) 


£  =  0 


BA  =  27.97° 


CB  =  83.40' 


The  resultant  mechanism  guarantees  the  satisfaction  of 
the  initial  problem  specification,  but  several  additional 
quality  factors  may  limit  the  usefulness  of  this  mechanism, 
To  further  evaluate  the  quality  of  this  mechanism,  an  anal- 
ysis must  be  run.   It  was  found  that  the  crank  AB  cannot 
rotate  a  full  360°  so  that  a  motor  could  not  be  attached 
directly.   Generally,  a  designer  would  stop  here  and  se- 
lect a  new  crank  from  the  locus  satisfying  the  initial 
motion  specification.   In  this  case  we  will  pursue  the 
analysis  to  demonstrate  additional  points.   The  analysis 
is  summarized  in  Figure  VI-3  for  order  k  =  0,1,2.   Note 
in  this  case,  that  both  positions  £=0  and  £=3  are  satis- 
fied without  reassembling  the  1  inks  i n a di f f erent  config- 
uration.  This  is  not  always  true  and  must  be  guarded 
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•-  -Y 


Figure    VI-2 
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0^  (deg) 


6i  (d«g) 


Figure    VI-3 
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against.   Note  also  that  the  position  specifications  for 
e^  are  not  satisfied.   Only  the  relative  displacement  be- 
tween the  £=0,3  positions  are  met.   In  other  words,  the 
same  locus  would  have  been  generated  if  Bgg  =  0°  and  633  ■■ 
18°. 


TR  Suspension 

The  state  of  the  art  of  automobile  suspension  has 
evolved  over  the  past  50  years.   Race  cars  and  passenger 
cars  have  mixed  objectives  of  handling  and  comfort.   Most 
of  the  comfort  problems  can  be  satisfied  by  the  proper 
application  of  energy  devices  such  as  springs  and  dampers. 
The  handling  problem  is  also  dependent  upon  these  devices, 
but  it  is  more  heavily  dependent  upon  the  suspension  geome- 
try.  The  position  of  the  roll  center  and  the  reaction  of 
the  body  mass  to  braking  and  bumping  have  lead  to  a  wide 
variety  of  suspension  geometries.   If  this  problem  were 
to  be  solved  with  spherical  motion,  several  advantages 
could  be  anticipated. 

Advantages  of  Tesar-Riddle  (TR)  Suspension 


Virtually  zero  tread  width  change  during  large 
bump  motion. 

a.  Reduce  scrubbing  to  zero. 

b.  Reduces  potential  break-away  due  to  scrubbing 
during  bump. 

Virtually  stationary  roll  center, 
a.   Most  roll  centers  tend  to  raise  or  lower 
during  bump  action  or  during  high  g  turns. 
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b.  The  location  of  the  roll  center  is  fully 
the  designer's  choice. 

c.  Stationary  roll  centers  should  enhance 
predictability  and  stability  in  high  g 
turns  . 

Zero  acceleration  of  vehicle  during  bumps. 

a.  Most  suspensions  change  the  distance  be- 
tween the  mass  center  and  the  tire  center 
which  creates  a  fluctuating  vehicle  ac- 
cel eration . 

b.  The  TR  suspension  reduces  this  chugging 
action  to  zero. 

Anti-lift  action  during  braking. 

a.   Anti-lift  requires  that  the  brake  forces 

create  a  moment  on  the  car  to  keep  the 

rear  down  during  braking. 


Solution 


The  Mercedes-Benz  450SEL  Touring  Sedan  was  chosen 
as  a  reference  body  to  which  to  fit  the  TR  Suspension. 
It  may  be  an  example  of  an  automobile  with  a  good  balance 
of  comfort  and  handling.   Its  standard  suspension  has  a 
fully  independent,  diagonal  pivot  rear  axle.   Items  1  and 
2  above  may  be  satisfied  by  arbitrarily  choosing  the  center 
of  the  sphere  to  be  six  inches  off  the  ground  at  the  op- 
posite wheel  (Fig.  VI-4a).   Items  3  and  4  can  be  satisfied 
by  requiring  the  wheel  to  also  move  on  a  path  that  is  a 
constant  radius  p  from  the  car's  mass  center  (Fig.  VI-4b). 
The  problem  statement  is  now  constrained  to  spherical  mo- 
tion.  However,  the  sphere  has  a  60-inch  radius  and  the 
wheel  motion  is  in  a  relatively  small  area  of  the  sphere. 
Initial  motion  specification  can  begin  by  a  planar  approxi- 
mation.  Assume  that  we  are  given  the  fixed  pivot  0^  (same 

a 
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position  in  existing  Mercedes  system),  and  the  direction 
to  the  center  of  gravity.   For  the  planar  approximation 
assume  the  axle  to  be  infinite  in  length.   Figure  VI-5 
illustrates  this  planar  approximation.   At  the  curb  con- 
dition, it  is  desired  that  point  C  move  perpendicular  to 


the  direction  of  the  mass  center.   Further  choose  0  C  to 
be  very    long  so  that  the  corresponding  point  on  the  in- 
flection circle  is  given  by  J^.   Arbitrarily  choose  the 
pin  joint  A  to  be  near  the  perimeter  of  the  rim.   The  in- 
tersection of  OgA  and  O^^C  locates  the  instant  center  P. 
The  Euler-Savary  equation 


J, A  = 


Ml 
OaA 


(VI-1) 


allows  the  calculation  of  a  third  point  on  the  inflection 
circle.   The  circle  may  be  drawn  through  points  P.Jg,  and 
J^.   Now  any  point  B  may  be  chosen  to  locate  a  new  con- 
straint that  will  replace  0  C. 

A  line  through  points  P  and  B  intersects  the  in- 
flection circle  at  J.  .   Rearranging  the  Euler-Savary 
equation. 


V  =  37b 


(VI-2) 


allows  the  calculation  of  fixed  pivot  O^.   About  six  quick 
iterations  on  the  choice  of  points  A  and  B  resulted  in 
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the  final  configuration  of  the  four-bar  O^ABO. ,   This  pro- 
cedure proved  valuable  in  predicting  the  spherical  solu- 
tion.  The  transformation  between  homogeneous  coordinates 
and  cartesian  coordinates  is  given  by 


A(zi,Z2)  ^  A(X,Y,Z) 

2    2  -■'/2 
Z  =  +  (1  +  zf  +  z^) 


(Vl-3) 


X  =  z^  Z 


Y  =  Z2  Z 


We  can  now  express  the  four-bar  O^CAO^  in  cartesian  coor- 
dinates as 


Oc{X,Y,Z)  =  (.2838,  .9478,  .1451) 

C(X,Y,Z)  =  (.1033,  0,  .9946) 

A(X,Y,Z)  =  (.1941,  -.0435,  .98) 

0^(X,Y,Z)  =  (.07,  .3333,  .9402) 

a 

Dyad  analysis  of  the  four-bar  will  yield  the  motion  prop- 
erties of  point  C.   This  motion  can  be  related  to  the  ro- 
tational parameter  a   and  g  by  Eq.  (III-3).   Further  analy- 
sis  yields  u  and  w.   The  corresponding  angular  velocity 
and  acceleration  about  point  C  are  given  by 


Y=(uj*c)-dsin6 


Y  =  ( 


u 


c) 


a    sin  3  -  ag  cos  6 


(VI-4) 
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The  time  state  motion  specification  is  now  complete  for 
the  3  ISP  case  (PPP).   The  curvature  transformation  yields 

B(X,Y,Z)   =  (.01324,  .0475,  .9988) 

0^{X,Y,Z)  =  (.2551,  .1769,  .9506); 

Fig.  VI-6  illustrates  this  solution.  It  is  felt  that  the 
TR  suspension  is  a  superior  suspension  in  accordance  with 
the  initial  objective.  The  most  significant  point  of  the 
solution  is  that  the  arbitrary  decisions  could  be  modified 
to  suit  an  experienced  suspension  designer.  This  example 
illustrates  the  creative  capabilities  of  wel 1 -f ormul ated 
analytics  translated  into  efficient  computer  programs. 


Poultry  Transfer  Device 

Like  most  industries,  the  poultry  industry  faces 
the  need  to  automate  processes  that  have  become  undesirable 
to  human  labor.   The  poultry  preparation  process  is  divided 
by  a  break  between  the  picking  line  and  the  eviscerating 
line.   Current  practice  connects  the  two  lines  by  a  con- 
veyor carrying  birds  with  random  orientation.   Workers 
pick  up  each  bird  and  place  it  in  a  shackle.   This  task 
is  messy  and  unrewarding.   Early  automation  attempts  have 
focused  on  picking  up  the  irregular  objective  and  rehanging 
the  bird  by  machine  hands.   The  uncontrolled  variation  of 
too  many  parameters  cancelled  the  reliability  of  such  an 
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attempt.   It  is  apparent  to  this  author  that  control  of 
the  bird  must  be  maintained  during  the  transfer.   This 
means  that  the  device  must  match  velocities  before  it  is 
released  from  the  picking  line,  and  transfer  with  matched 
velocities  of  the  second  line.   The  object  of  this  design 
will  be  to  synthesize  a  mechanism  to  transfer  from  the 
picking  line  to  a  transfer  slide.   The  same  mechanism  can 
be  used  on  the  other  end  of  the  slide  to  transfer  to  the 
eviscerating  line. 

The  hock  cutter  is  the  last  operation  of  the  pick- 
ing line  and  produces  1.8  birds/sec  with  a  linear  velocity 
of  10.8  in/sec.   Inspection  of  the  machine  shows  that  a 
40-inch  sphere  would  keep  the  volume  in  proportion  to  the 
frame  and  maintain  an  approximate  straight  line  while 
matching  velocities.   The  bird  will  be  delivered  to  a 
slide  by  an  in  and  out  cuspidal  motion  (Fig.  VI-7). 

The  5  MSP  motion  specification  set  to  describe 
the  desired  path  is  given  by 


point  c 


point  a 


I 

k 

a 

e 

Y 
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2 
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0 
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0 

0 
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1 

0 

0 

.54 

where  the  first  three  positions  describe  the  cusp  and  the 
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transfer 
slide 


Figure    VI-7 
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last  two  describe  the  matched  velocity.   The  curvature 

transformation  yields  six  real  solutions  of  which  any  two 

may  be  chosen  to  form  the  four-bar  constraint.   The  chosen 
constraints  are 


A(X,Y,Z)  =  (-.9046,  -.2563,  .3406) 


B(X,Y.Z)  =  (-.6259,  -.6441,  .4398) 


C(X,Y,Z)  =  (-.7070,  .7072,  0) 


D(X,Y,Z)  =  (0,0,1) 


AB  =  28.23' 


BC  =  90.74° 


CD  =  90° 


Since  CD  =  90°,  this  mechanism  is  the  equivalent  of  the 
planar  slider-crank  mechanism  (Fig.  VI-8).   This  device 
will  constrain  the  specified  point  path  but  it  is  not 
yet  coordinated  to  the  input  drive  of  the  picking  line. 
Analyzing  the  four-bar  ABCD  will  yield  the  position, 
slope,  and  curvature  of  any  point  b  along  the  path.   Ar- 
bitrarily choose  a  point  b  along  the  path  from  the  pickup 
to  release  and  choose  a  relatively  slow  velocity  magnitude 
in  the  direction  of  the  slope  at  that  point.   The  cuspidal 
motion  is  satisfied  by  the  four-bar  geometry,  so  the 
second  velocity  point  can  replace  the  cuspidal  motion  of 
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point  c  as  a  design  position.   Analysis  will  yield  the 
angular  motion  of  link  AB .   This  motion  is  to  be  coordi- 
nated with  an  input  crank.   The  problem  is  now  reduced 
to  a  function  generation  problem  where  the  solution  must 
yield  a  four-bar  drag  link  so  that  a  continuous  motor 
drive  will  provide  the  input.   The  cycle  time  of  1.8 
bird  per  second  translates  into  an  angular  velocity  of 
the  input  crank  of  11.3  rad/sec.   This  problem  can  be 
solved  by  any  of  the  mechanisms  listed  in  Table  4  and  by 
any  procedure  in  Section  4  of  Chapter  III.   We  will  seek 
a  drag  link  by  choosing  a  short  base  link  (AB)  of  ten 
degrees.   The  angular  parameters  of  the  input  link  (GH) 
and  the  analyzed  link  (AB)  complete  the  necessary  data 
for  problem  specification. 
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units 


0 

0 

3.14159 

1 .0724 

rad 

1 

1 

11 .2923 

-2.344 

rad/sec 

2 

0 

4.1888 

0 

rad 

3 

1 

11.2923 

-2.344 

rad/sec 

This  data  set  can  be  used  in  the  APL  function  LINK3  to 
yield  the  locus  of  link  6F.   The  coordinates  of  the 
four-bar  drag  link  solution  are 
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H'(X,Y,Z)  =  (0,0,1) 

A'(X,Y,Z)  =  (-.0868,  .1504,  .9848) 

F'(X,Y,Z)  =  (.0870,  -.0256,  .9959) 

G'(X,Y,Z)  =  (-.2244,  -.0898,  .9703) 

The  relative  location  of  the  references  is  arbitrary  for 
function  generation.   This  allows  the  transfer  of  the  de- 
vice (Fig.  VI-9)  so  that  points  A  and  A'  coincide,  and 
points  B  and  F  are  in  the  same  moving  shell.   Hence  the 
dyad  HGF  drives  the  four-bar  ABCD  with  the  proper  dynamic 
motion.   The  solution  shown  in  Fig.  VI-9  could  be  further 
optimized  by  repeated  computations  of  the  above  procedures 

Conclusi  on 


This  work  has  met  the  major  proposed  objectives. 
The  analytical  tools  of  spherical  synthesis  and  analysis 
have  advanced  to  a  level  of  understanding  approaching 
that  of  planar  kinematics.  The  past  works  of  Drs.  Tesar 
and  Duffy  have  allowed  a  coordinated  mesh  of  planar  syn- 
thesis and  spatial  analysis  to  yield  comprehensive  tools 
of  spherical  mechanism  design.  The  combined  works  of  D. 
Tesar,  J.  Duffy,  H.  Dowler,  and  D.  Riddle  have  accomplished 
the  following  objectives  of  spherical  kinematics: 
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Figure    VI-9 
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1.  Curvature  transformation  for  3,  4  and  5  mul- 
tiply separated  positions 

2.  Time  based  curvature  transformation 

3.  Colaminar  Inversion 

4.  Classification  of  twelve  solvable  problems 

5.  Synthesis  Procedures  of  four-,  five-,  and  six- 
bar  mechanisms 

6.  Closed  form  angular  displacement  analysis  of 
all  four-,  five-,  and  six-bar  mechanisms 

7.  Time  state  analysis  of  translational  and  ro- 
tational motion  up  to  the  k-th  order  for  all 
four-,  five-,  and  six-bar  mechanisms 

The  design  of  spherical  mechanisms  will  reach  a 
parallel  level  with  planar  mechanisms  after  the  completion 
of  two  additions  to  this  work.   First,  a  comprehensive 
user's  manual  of  interactive  APL  functions  is  needed.   A 
few  of  the  functions  listed  in  the  Appendix  could  be  fur- 
ther generalized  (e.g.  LINK3  is  coded  only  for  the  case 
PPP-P).   Most  of  the  existing  functions  make  efficient  use 
of  the  central  processing  unit,  but  greater  care  could  be 
observed  in  storage  requirements.   Second,  the  analysis 
tools  for  systems  with  mass,  springs,  and  dashpots  could 
be  added  for  general  dynamic  analysis.   This  would  lead  to 
a  basis  of  shake  balancing. 

Optimization  is  the  most  significant  missing  ele- 
ment in  spherical  design.   It  is  also  needed  in  all  other 
areas  of  machine  design  including  planar  kinematics.   Planar 
optimization  is  currently  being  examined  by  several  re- 
searchers.  Classical  methods  generally  require  the 
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formulation  of  a  multi- parameter  cost  function.  This 
author  believes  that  this  is  unrealistic  due  to  the  num- 
ber of  dependent  variables.  An  alternative  approach  se- 
quentially filters  undesirable  mechanisms.  The  sequence 
is  ordered  according  to  computational  efficiency.  It  is 
suggested  that  this  philosophy  will  also  apply  to  spherical 
design. 


APPENDIX 


INTERACTIVE  SYNTHESIS  AND  ANALYSIS 
COMPUTER  PROGRAMS 
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7    Zff/  ABI  A 
[1]         Z-H(p/l)pO 

[2]         ZCl    H    ;23-^"lOi4Cl    4    ;1] 
[33         ZCl    4    ;l]-{-/lCl    U    ;2]^20ZC1    U    ;2])    VATAfJ   All 

H    3]*20ZCl    U     ;2] 
[U]         Z[2;2]^4[2;1]t0o/1C1;1] 

[5]         ZC2;1]^(/1C1:    2    3]-.Xi4[2;    3    2])44[1;3]*2 
C6]         ZC3;2>((i[3;l]x0Oyl[l;l])t4[l;l]x;lC2;i:*2)vl-;lCl;l3* 

C7]         ZC3:l]^((/l[l;3]x4[lj     2    3]-.x;i[3;    3    2])-2x4C2; 
3]x;l[l;    2    3]-.x-4[2;    3    2])MCl;3]*3 


V    Z-eB    AUG    Cfl\A\X\AV 
[2]         ^-iVpO.BC;3].(-B[;2]).BC;    4    5    fi]  .CC  ;  2]  . -C[  ;  1  ]  .  C 

C3]        i:4+(JpB[;    4    5    6:.^r.;    W^ll^l^.^'fl]\l\\   gr. 
1.1]    BC;l])  +  /VpC[;    4    5    6],0,C[;3].(-C[;2]),flL, 

2i,-S[jl3.Cl-l]0  .       nr 

[U]         ;f-H(7^-    0    0    2)p(0.0.[1.5]    sr;5])  +  (0,CC;4].C 
1.5]    B[;6])+Z?[;U].0,[1.5]    ^C;5] 

[5]         Z^X    QUADIV   A  .  x     /       ^r    m    r 

[6]         X^(N-    0    0    2)p(0.0.[1.5]    3[;8])+(0,CC;7].[ 

1.5]    Z?[;9])+iJ[;7].0,[1.5]    CLi9l  ^r     ,  t    r 

[7]        4Z)^A7pO.Z?C;6],(-BL;5]).BC;    7   8    9].^C;5].-(7C;4],C 

[8]         AhtAD^iNpBZi    7    8    9].CC;    7    8    ^^  •  ^  "f  ^^  ;f  ^  ^  1°  '  ^^r 
1.1]   flC;4])  +  f/pC[;    7    8    9].0,i7C;6],(-^[;5]).BC; 

5],-BC{H].Cl.l]    0  .     -,     ,    r 

'  C9]         Z-HZ,(;f-(A'-    0    0    2)p  +  /i4PxZ,r2]    Z.L 

1.1]     Z)    Qi/>1DIV   A  n  „s 

[10]       Z-Z.((2[;i3]+.*2)*0.5).C1.1](Z[;    U    5    6]+.*2)* 

0.5 


V  B"rA    CRANK    XiAB  iW  ;T  ;N  i8',E 
C2]    T^(RloilQ0)^Xi32+Xi5l^'l-i-\N'>-l  +  l(Xi^l-Xl3l)^XL 

5] 
[3]    B'^iAB<rRxXZll),V'^R''XZ2'] 
[H]    g^((10>lZ?)xl02').(-(10/li?)x202'),C 

[5]    3^1](tI^I    3  2  «?(.V.  3  3)p2'.(-x/^Cl  2].*r).(-xM[ 

1  3].*r),o.(^C3]Tr).(-/iC2]*r^ooyi[i]).3+;i)x2.[ 

2]  5_,Cl.l]  2  .    .^    r  .-,  r 

[6]   fl-HS,+/rx,'?.[2]  ;?.[!. 1]  /?^,Vx(-SC;2]).2C;l].L 

[7]    B^S.+/2'xi?,[2]  .7. [1.1]  .^-Hf/xC-/?:  ;2]),.7[  ;1].C 
1.1]  0 
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V    C^A    DYAD    X',B;DiTiP;BCBD;CD:R;KiN',Mi'!2',:!^;Cl 
[1]         K^(l'C->-P^X[ll)  ,iCD*-Xl2l^R^Oil^0)  ,^'-'\>    3    3     ,'^^pDD^ 

~20  +  /(,B^Al;    1     2    3])xP-e^[;     10    11     1?] 
[2]         -*RED^\OeK^l^\T^i  (  2oCD)  -  (20BC)y^20!.D)  ^  (loBC  )<loBD 
[3]         pt^\(Qi2)-'20((,20BD)-(20BC)^20CD)i-il0BC)^10Cn 
[4]         P-e+/C2](Z).i7,[2.5]    D)^    1    3     2    <5)A/3^;/p/? ,  ( -F[  ;  1  ]  x/?r  ; 

2]4i?),(-S[;l]x;7[;3]iiO,0.(:/L:3]T.?).(-Sr;2]^i7^0OZ?[;l]) 

[5]         .'.:f4W/p(-P[  ;2]i/?),(?r  ;l]v/O.0,(-P[  ;     1    23t/?,[ 
1.1]    i?-»-(PC;     1    2]+.*2)*0.5).0,0,0,r 

[6]         C-Ht/[2](/?.P.C2.2]    Jl^(ilO-2oT)^loBC)  A-T^IODC)  ,1 

1.1]    20i?C)x'./^    PDT    /^3 
[7]         C1'^-  +  /C2](C.^,C2.2]    r7)-(l     3    2    ^NpAL  ;    4    5    G    13 

14   15], 0.0. [1.1]   /^po)   .'T  1   3    2  fi^r-H.'^i.vy  ;.'pii.I?. 

2xC 

[8]  C^C,Cl.-+/irxP.C2]    .'?.[1.5]     ?^(+/(.Vp.tC;     7    8     9 

16    17    183.0.0.C1.1]    ."pO)xC",[2]    C.L 

1.1]    C)  +  2x  +  /(/;p/l[  ;     4    5    6    13    14    1  5  ]  ,C1  )xCl ,  [  2  ]    Cl.L 
1.1]    CI 

[10]    PFZ):H.i.(.yl^(CZir?5E'-*-;O/>l).0tn^'    ***    CLOSURE   P.FDUCTION    ** 
*' 
V 


7  Z-t-MINV   BiMil 

[1]  Z-t-    0    3    3    p<^/-*-(  J-^l)tp5 

[2]  Z^Z,[1]EBCI;;J 

[3]  -►2xA/c:i"-i-r+l 

V 

7    Z+-M   PD?    ;•/ 
[1]  Z-^((2+p;*/)p     0    "1    "2)c|>(+/."xA').(+/"/xlbC2]    N),L 

2.2]  +  /./x-i©[2]    A'^    13    2    fi)A' 
7 

7  R-A    QUADIV   B;D\EiI 

Cl]  D-^oByJ'  l  ,OpP^(0,(p/  )C2])pO 

[2]  P-Hff,[ij,((~2+p/l)pyl[I;;])K('2+Z))p5[I;;] 

[3]  -*2xiOC  1]>I-^I  +  1 

V 

7    Z-'rRELA    A\B;I 
Cl]         S-*-"20(+//lC  ;9  +  I]x;l[  ;18+J]),(+/4C  ;7]xylC  ;18+J]),C 

l.l]  +  //l[;J]x/lC  :9+I-Hi3] 
[2]  Z^"20(-(205r;2])-x/20ii[  ;     1     3])^x/lo3[;     1     3] 

[3]  Z^(02)  |Zx"l*0>  +  //ir  ;9+I]x( -//.[  ;     20    2l]x;ir;     3 

2]). (--//}[  ;     19    2l]x/lC  ;     3    1]),C 

1.1]-/4C  ;     10    20]x/l[  ;     2    1] 

[4]        S-»-*x/io5r;    1    3] 

C5]  Z-»-Z.[l.l](5vl02)x  (+/>![  ;3+J]x,1[  ;ia+J])  +  +  //lC  ;I]x/l[  ; 
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[5]         Z-'-Z,[l.l](BTlo.")x(t//l[  ;3tl]x.4[  ;  18+ J]  )  +  +  //.[  ;I]></1[  ; 

21+J] 

[6]  Z^Z,(-(2C  ;2]*2)^30Z[  ;1])  +  (Z?M0Z[  ;!])'<(  t/4[  ;6  +  J]x/!C  ; 

18+J])  +  (2x  +  />5[  ;3+J]x/?C  ;21tJ])^  +  /ylC  ;i']x4[  ;2U+J] 
7 


V    Z^STItlV  A;MiB;C;K;Li:n\M^',V',R 
R-^-A  [  1  ;    1    2  ]  +  .  *  2 


[  1  ]         R-^-A  [  1  ;    1    2  ]  +  .  *  2 

[2]         M^{B  ,C,OA-C'>-iAll;ll^Al  ;  2  ]  ) -/IC  1  ;  2  ]x/l[  ;  1  ]  )  ,  (B-H(y?[  l  ;  l  ]x 

/IC  ;  1  ]  )  +4  [  1  2  ]  x/1  [  ;  2  ]  )  .  0  .  0  .  0  ,  [  1 . 1  ]  (  ;>1  +P  -4  )p/?  )  *i? 
[3]  W^(X-h(6xA'),     3    3)p     2    1     3     4    fij(6,j^:,     3    3)pM 

i'42  Z^+/Wx     13    2    <S?;fp2?,Z),C2.5]    P-'-AiiL-^    1     2    3    10    11 

12    19    20    21    28    29    30    37    3R    39    46    47    U8] 
[5]         //3-H(B,C.O,(-(7^(>iri;l]x/lC  :5])-i4[l;23x/ir  ;  4]  )  ,  (B-.-Cyir  1  ;  1  ] 

X  >1  C  ;  4  ]  ) + yl  [  1  ;  2  ]  X  yl  [  ;  5  1  )  .  0  ,  0  .  0  .  [  1 .  1  ]    i:  P  0  )  i /? 
[6]         .'<f3-*-(^^(6x^),     3    3)p     2    1     3    4    ^{S,K,     3    3)pi"f3 
[7]  Z-<-Z.{+/W3x£))  +  +/.:'xl/.K    1     3     2\KpV,V,l 

2.2]    V-<-/l[;L  +  3] 

W4*(S,C.0.(-C*-(/li:i;l]x/}[  ;8])-;lCl;2]x;}C  ;7j).(B-e(/lCi;l] 
x^C  ;7])+4[l;2]x/ir.  ;8}),0,0,0.[1.1]    ;rpO)*i? 
MH'*-(K'>-(6^K),     3    3)p     2    1     3     4    6?(6.A',     3    3)p?/4 
Z-*-Z,(+/'/4xZi)  +  (  2x  +  /-'3x7)  +  +  /.7x    13    2    <S(KpA  ,A  ,1 

2.2]    /1^>1[  ;i^  +  6] 


[8] 

[9] 
[10] 


[11]       Z^    2    1    3    <Si(K,     6    9)pZ 


7  Z^STl    AiX5',YSiZSiX2iY3;Z3;X51iYSliZ51;T 

[I]  A:5'»-1x.o/1[4    8] 

[2]  y5-f--(l     2    x.>l[3    4])+    2    12    x,0/l[3    4    8] 

[3]  Z5-^(2x.Oy4r3    4])-     1     1     2    x.o/l[3    4    8] 

[4]  ;f3-^-lx.oyl[6    7] 

[5]  y3^-(l    2    X.O/1C2    6])+    2    12    x.o4[2    6    7] 

[6]  Z3-t-(2x.o>}C2    6])-    1    1     2    x.o<4[2    6    7] 

[7]  2'-»-(o*lfl)x~l  +  i36 

[8]  ;!r51-^(-^5x202')-y5xioT 

[9]  Y51'-(.{2oAll2  )^Z*-(XS^10T)  +  Y5^20T)-Zbxl0ALlZ 

[10]  ZBl-nCdO^Cl])'   ■)  +  Z5x20/l[l  ] 

[II]  D-^(  ;sr3  x;f  5  1 )  +  (  2C..U  5  J  )  -  (  y3  X  i'5  1  )  +  Z3  X  Z5  1 
[12]  B-H2x(;!r3xy5l)+^51xy3 

[13]  ;^^(y3xy5l)  +  (20/J[5])-{;f3x75l)  +  Z3xZ51 

[14]  K'^(X5-^1E~2S<K)/  -^(3*2  ) -^xD^K 

[15]  Z^(.X5/T),(^(2.r/.:5)-^2-'3o{(B,B^X5/-B)i-K,-K'*-K* 
0,5)i2^A,Ai-XS/D 


7    Z-t-A    XPR    B 
[1]         Z-e((;l[2]x5[3])-/l[3]xZ?C2]).(-(/iri]x5[3])-;l[3]xB[l]).(^ 

[l]xS[2])->l[2]x??ri] 
7 


174 


V    //  ASSIGN    L',A  ;B 
[1]         ->(3    i|    =N)/L1,L2  ^     , 

-hA/[ L-^  ;  2  ]  )  .  (  CA^2oA  )  ,5 A^\oA<- !■".{.  L-!l  ;1'\ 
[3]         -.0xil  =  /y.0t(41^';[L+l-^;l]v/l).-^l*--'CLtl-A^;2]*/l^.VCLtl-A'; 

3] 

[H]  ->0.(S2<-(/<f[3;2]-Blx/l)ifl),/12^(A/C3;l]-.ll>«.1^.'/[3; 

3])T5-e//C2;3]*2  ^  ^     ^ 

[5]       LI  :-*0  .  (SS-^C  /■/[  4  ;  2  ]-  f /}  Ix.'/C  H  ;  3  ]  )  +  3x1?  2xF  )  t.1  )  ,4  3^(//L 

'+;l]-(/llx;.f[i+;3])  +  3x/12x3-«-x/;fr2    3    ;  3  ]  )  */l^A/[  2  ;  3  ]  * 

3 

[6]       L2  :  -vO  .  ( /!  4-.'''C  5  ;  1  3  )  ,H  4^'f  C  5  ;  2  ] 


7    Z-<-C>}5F    ^;i'/ 
[1]         K-^    1    2[l+'P''=A'].Z^iO 
[2]         Z-^2,~l  +  i~l+A'-*-i'x2 
[3]         -^2xO<p^-^ff+A' 

V 

V    Z-^CUDIC   XiAiBiCiD 
[1]  -».Hxi(Z  =  0)vO^A'^((2^(/lxC)-??*2)*3)  +  (D^(-ZJ*3)  + 

0.5x(3x4x(S^A[2]v3)xC^;f[3]T3)-/lx(4^1+j:)xD-e    l+;f)* 

2 
[2]  Z-h(3    1    p((2xxD)x(Z*0.5)x2O((  ^3)x    20(  [  Z5  )  ^  (  Z-f-Z  )  * 

1.5)+  0  1  2  xO243)),0 
[3]  -►O,  ,Z-t-(Z-  3  2  pS,0)t/1 
[4]  Z*-    3     2    p(T>l)x(-ii-  +  /^).0,Z,Z),(Z-e(-(42)x  +  /C)-S).-O^(  + 

2)x(3*0.5)x-/C-»-(Z?+;f,-Z-<-I*0.  5  )**3 


7    y-t-fj;ir    ;ir;WF;S;l;5B;C>J  ;CJ?;S(7;C(; 

[1]         0    ASSIGN    1 

C2]  MF^    3    3    p(C3xCC),(-i7Z?x5^),5/?.((^/-x5^)+5/!x55xCC).((C>lx 

CO -S/l  X53X5/7  )  .  ( -iM  xff;? )  ,  (  (  5yl  xSC  ) -C>1  x5S  xC(?  )  ,  (  (  54  xCff) +C4 

nS3i(SG)  ,CA*C3 
[3]  y-^fi?//F+.xZ 

7 

V    A-^-G'fL    KiAl  iA2  ;A:i  iAmUl  iB2iB3  ',Bi*  iSA  iCA  iS5  iCB  iSGiCGiS  iT 
Cl]  0    ASSIGN    L^{K=0)/\pK 

[2]         4'^((p/:)»9)p0 
[3]         ALL-.l^il-CB^CA)  AiSGxCAy(SB)fSA^CG)A(CAxSB^CG)-Sa^SA) 

,(~l+Ci3xCC)  ,(C2;x55),(J/lxC".'?),5B.((-5Bx5.4x5r;)-l-C/lxCff), 

[l.5](C4x5(7)  +  r^x5'!?x5/l 
[4]         4CL;]'H4CL:]-((pL),0)p//0-^4[l:] 
[5]         ->(9y?'xi~l€^ 
[6]  1    ASSIGN    L^{K=l)/\pK 

[7]         4CL;    1    2]-<-((41x54xCS)+SlxC4x5."),C 

l.l2(SGx(.SA■^(Al  xJA^Sn  ) -Bly<CA^CB  )  )  +CGxCA^SB  +  Al 
[8]         AIL;    3    i4]-e((S^x(;4x-5/?+41)+CCx(-54xl+41x5S)+/?lxC4xCi?), 

[l.5]-(5ffxC5)+C^xS2?x51 
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[9]         ALL;    5    6    T]^(  (CByCO-SB^SG^Bl)  ,  i  i -^^'^''2A^SB)-^A1^CA^CB  ) 

r  1  1 1  Bi  X  '^^ 

[10]       all';     8    q]^((S,'7x-C)  +  C^x/;),Cl.l](5Gxn^-5/lx5J?+/}l)  +  Cf7xC^( 

Bl  x5^  X  CB  ) +C/1 +yl  1  X ^/ xS2^ 

[11]       ■*0UT^\~2€K 

[12]       2    ASSIGN    3  v  ^    / 

[13]       >l[3;l]^(5/lx(-2x41xBlx5B)t/12xCB)  +  (B2xC/lx57?)  +  Ci4xCBxC'»-( 

[14]       4[3;2]-H(5(;x-(/l2xiMx5B)  +  (2x/?lx/?lx5/lxC5)+(CXx5F)  +  ((7xC/lx 

5B)  +  (2xyllxCi4)-B2xC;lxCB) 
[15]       ;l[3;2]-K^[3;2]+Cffx-(2x/llx5/lx5B)+6',1  +  (5^x^l*2)-( 

2xBlx<:7/lx(:g  )+A2xCA 
[16]       /l[3;3]*(SGx(5>lxl+yll*2)  +  (2x/llxS/lx5B)-(42xC^)  + 

2xBlxC/!xCB)+CGx(B2xC/lxi7B)-(2x>}lxC^)  +  (yl2xS/lx5B)t( 

2xy41xSlxSi4xCB)  +  C^x5Bxi  +  C 
[17]       /1[3;     4    5     6    7]^((5(;x2xBlx.7B)-C(;x(CFxl+Bl*2)+S2x5B).(( 

5Cx-(B2x5B)  +  CBxl  +  Bl*2)-2xfllx5BxC(;).(-(5/1x(B2x5B)  +  CxC5 

)+C/1x(2x^1xS1x5B)-;12x(:B).(-5BxS1*2)+B2xC.'? 
[18]       ^[3;     8    9]^((5ffx-i:)  +  C(;xr),(S^x2'-^-(2x51x5/lxf7B)  +  (yl2x5/l)  + 

(2x/llxC>lx5B)+!7;txl+/l*2)+C*Gx.7^(-J4x5Pxl  +  C)  +  (S2x5>lx(7B)  + 

( -2x/}  1x5/1  )tr/lx(>12x5B  )  +  2x/?lxBlxCB 
[19]    OUT-.A*-    1    0    +;1 
V 

7    Z^K   MSP   MiAiN  i3  ;C  iD  iE  iF  iG  iH  iA  ',P  iQiR  iQ  irlO  iA!fiI  iJ',L  i?  ',U 

;21;Z2 

[I]  AM^A^^TML    Gl'lL    A'.^-^-i  QxJ-hI 

[2]      Ll:fl^(A   >1[3    9    i/VD.CA   >1C8    2     ;//]). (A    >1  [  9    2     ;/l/])  +  A    /1[ 

38     ;A?-t-    2112323    3[J,J+4]] 
[3]         S-»-fl,((A   >1[1    9    ;^])  +  A    yl[5    3     ;;V]),((A   Ai2    6    ;;V])  +  A    Al 

1  8    ;A^]).A    4CG    1     ;/n 

[4]         C-^CA  i4[3    4    ;iV]).(A    /1[2    5     ;A']),(A   4[4    2     ;W])  +  A    /1[ 

5    3    iN-\ 
[5]  C-H^,((A   /1[1    4     ;/V])  +  A    /l[3    7     ;//]). ((A   AiS    1     ;fl?])  +  A    Xl[ 

7    2    ;JV]),A   yl[l    7     ;fl?] 
[6]         D-^iL   4[4    9    ;//]), (A   /1[3    5     ;//]). (A    >1[9    5     ;yV])  +  A   4[ 

4    8    ;iV] 
[7]         D*-DA(i^   AL6    4    ;//])  +  A    All    9     ;;V]).((A   /1[5    S    ;W])  +  A    /4[ 

7    8    ;iV]),A    Alls    7     ;//] 
[8]  ^♦•+/(6     3    p/l[4    9    3     ;!■<-    14    4    4[  J]  ]  )  xB  , -C ,  [ 

1.5]    D 
[9]         H-^//.[0.5]+/(6    3    pAl2    5    8     ;     ])xD,-B,[ 

1.5]    C 
[10]       //-.-//,[  l]  +  /(6    3    p.'l[7    5    1     ;I])xS,C.C 

1.5]-0 

[II]  P-^B[2;2],(//[2;5]  +  ;/[3;2]),(ffC2;6]+//C3;5]),f?[3; 

6] 
[12]      d^d . //C 1 ;  1  ]  .  ( ^^t  1 ;  3  ]  + f/r  2  ;  1  ]  )  .  ( ff  [  1  :  4]  +^[  3  ;  1  ]  ) .  ( //C  1 ; 

2  ]  + // [  2  ;  3  ]  )  .  ( // [  1  ;  5  ]  + ;/ [  2  ;  4  ]  + // C  3  ;  3  ]  )  ,  {  ff  [  1  ;  6  ]  +  B  [ 

3 ; 4]  )  t? 
[13]       -*-l(3=pK)  ,5>JW  +  l)/F0Un,Ll 

[14]       A^^    4    10    pd 
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[15]       e7-^-43p23'4l34124123    ^I^l+F^O 
[16]    L2://-^4>(A   /!['+    1    2     ;P]),((A   >1[4    1    3     ;P])t&    A[5 

1    2    ;P]).((A   /1C5    1    3     ;P])  +  A    /1[6    1    2    ;P]).A    /![ 

6  13     ;P-^J[I;]] 

[17]       F--£'+(""l*I)'<(4[7;I]x~l+/y),(4C7    8     ;I]+.x    2+//)»(/l[ 

7  8    9     ;J]+.x"3+W),(/l[7    8    9    10     ;I  ]+ .  ^'.'^  )  ,  (^  C  8 

9    10     ;I]+.x3+ff),(/li:9    10     ;J]  +  .x2+//).>lClO;I]xlt;/ 

[18]       -•■L2XI  5>I-^J+1 

[19]      ->0xZl--il>fl/^p22^(0  =  Z2[  ;2])/(22-^P^/?7'    F)[;l] 
[20]       Z-*-((W,3)p4[l;J]).[I-»-l](2l/lC2    3    ;J]  )  ,[  1  ]  (  Zi/1[ 
4    5    6     ;J]),C0.7](Z^(//,l)pZ2)i/l[7    8    9    10     iJ^x 

[21]       Zl-^-Zl.-CA    ZCl    2    4    ;I;])*A    ZCi3;J;] 

[22]       -»-("l+l2  6)xi;i/iI^I+l 

[23]    L3:Z-.-(ZxZl),[l](ZxZ2).[0.5]    .-^^C  1  + (  Z  1*2  )+Z2*2  )  * 

"0.5,0t,Q-'-    0    2    pJ-f-    0    0 
[24]      Z?-H(--/zr  ;J]x/l.-fr3    2    1     ;J]).(-7'Z[1    3    2     ;^]x^'/[ 

9    6    8     ;J])»Cl.5]-/Z[  ;J-J  +  l]x/l/f[4    5    7     ;I-^    2    3] 
[25]      H.(~i  +  i26)xi//>l+p^^r).[l]4>(-DC  ;1])!?PC  ;    2    3] 
[26]       ->O..Z^(FJI    2),(?'^(^xZ?,[1.5]    P).Z}^(l+(?+.*2)* 

~0.  5 
[27]    FOUn:X^(-/A!:LS    7     ;l]x4'/[G    8     ;  1  ]  )  f /)-h  +  M,'<[  u    5     :l]x,4:'/[ 

8  9    ;1] 

[28]       y-<-(  +  /4;/[7    4    ;l]x4.\'r9    f>     ;l])*i5 

[2  9]       C-H(4Cl]xZ?*3)  +  (il[2]xZ;x/)xF)  +  (^[4]xExFxD^2or)t<i[ 

7]x£;xEx£'-»-lor-<-(Of  1  fl  )x  il  R 
[30]       C^C,[l.l](dC3lxPxD)  +  (dC5]xDxJ.:)  +  (/l[8]x/7x.»?)  +  (4[l]x 

3x;fxDxZ))t(AC2]x(2xA'xDxF)  +  JxPxZ))  +  (ji[i+]x(2xyxZ?xP)  +  A'xFx;r 

)  +  3X(d[7]xyxFx;? 
[31]       P-H(4[e]xZ))  +  (yl[9]x,'?)  +  2x4[8]xyxF 
[32]      C-^C,.'?  +  (3x4[l]x;ifx;fx,9)  +  (aC2]xJx(^xF)  +  2xyx.r))  +  (4[ 

4  ]  X  yx  (  y xl? )  +  2 X IxF  )  +  (  3  x,J. [  7 ]  X  yx yx '^  )  +  (  2  x4.[  3  ]  xZx /7  )  +  ( ^ [ 

5]x(JxF)  +  yxZ)) 
[33]       R^(0.5iC,C*-I/Cl',ll)^{D,D^-I/Cl  ;  2  ]  )  tZ?  . -0^(  (  J«hD20  )/Z}-h(  <7 
[ ;2]*2)-4xx/C[ ;     1    3])*0.5 

[34]     2i"»-(i?x2or^r,r-Hj/2')+;f 
[35]     -^z,3,;;-^pZ2-H(ifxio?)  +  y 

V 

7    Z^PBRT    PiBidN  iTOL'.AiD'.I 

[I]  rOL-^-lf^lO.Ot.Z-e    1    2    pO 

[2]       ■♦(5   4  3  =pP)/rT,rri,o.D 

[3]         Bf-C^iN-^-pA^    0    0     ,-)pO 

[4]  J-H3,0+i/^(l    "l     x"2ti4)4lt"3+/l 

[5]  -►5xi//>J-Hl+l,0tCCI]^(UCJ]-»-^Ci']  +  ^'/+.><-^[-Z'-    1     2])  +  f/+.xC[I- 

1    2] 
[6]  U-^Ui-D-^iiBLN-    0    l]-.xC[;/-     3    2]), /?[//-    1    0]-.xC[/;- 

1   2])i(:cw-   2  i]-.x^[;;-   2  3] 

[7]         -^(2+J^3)xxr{?L<r/ 1^ 

[8]  Z-<-Z,[l]    nUAD    U^l,-U 

[9]         -►2.P^?    PDIV    U 

[10]  f?D:-^0.,2-»-    1    0    +Z,[1]    OlfAD    " 

[II]  C5:-^0,,Z-»-    1    0    +2,ri]    CUBIC    P 
[12]  Q2':-^0,,Z^    1    0    +2.[1]    QUA^IT   P 

7 
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V  Z->rQUAD    X 
[1]         X*XiXill 

[2]  -»-i4xiO>Z-*-(HxA'[3]  )-X[2]*2 

[3]         H.0.,Z^(2    1    p(f2)x(-Z[2]-Z*0.5),-.Y[2]  +  Z* 

0.5),0 
[4]  Z^    2    2    p(^2)xA',(-Z),(X^-lC2]),Z-H-(-Z)* 

0.5 
V 

7    Z^QUART    XiY',RxA\B',D',E 
Cl]  Y^,l^, CUBIC    l,(-X[2]),((x/;rCl    3])-4xZ[4]),(-x/.YC 

1  1     H])  +  (x/4,,y[2    4])-x/(A'"el4.ZiZCl]  )[3    3] 

[2]  -^4xi(i?A.=0),0+>l-f-(((Z[l]x;sf[l]x3T4)-2xAf[2]),0)-/?    C?    -'^'^C 

5>0)c{.0,(  |5^y-A'[2]-^Cl]xA[l]*4)* 

0.5 
[3]         -v5,Z)^(*2)x2  +  ,2    £/?   ^+/?-h((x/u,XC1    2  ]  )  -  (  8xZ[  3  ]  )  +  XC  1  ]  * 

3)    £Z5    nx/? 
[U]         Z)*-(^2)x2  +  ,2    £ifi4  +  i?-^2x(7>0)c{)y^(y*2)-UxAr[u])* 

0.5 
[5]         ffH-(i2)x2  +  ,2    £i?    /1-Z? 
[6]         Z-H    U    2    p(>l+J^+/7),(>^  +  .7-D),(yl-i?-ff),(/l^(-Ar[l]Ti+),0)-(i?-<-i?-f 

2)+F 
V 

V  GML^TML    GMiG;L 
[ 1 ]         GML^GM 

[2]         -^Oxi~l€;s:^l+A' 

[  3  ]         GMLiL  ;  ]^G//[  L;]x«}(9.pL)p(  C^l  +.'/[  ;  3  ]  )  [L*-(  .ir  =  1  )  /  i  p  A'] 

[U]  ->0xi~26/^.' 

[5]         GULL  2  ;  ]-»-( CJ/VC  2;l^G[.l  2*2)+G!fl  1 ;  ]  x£;[  2  ] 

[6]         -*-0xi~3e/: 

[7]         Gi'lLL  3  ;  ]-  (  C./C  1  ;  ]  x.'7[3  ]  )  +  (  3x^,'/[  2  ;  ]  x  x/(7[  i    2  ]  )tC/VC 

3; ]xCCl]*3 
[8]         -►0xi~4£;i: 
[9]         ffWLCif  :]-^(c;.'f[l;]x^[i+3)  +  (C7C2;]x(i+xx/(7Cl    3])  +  3xffC 

2  ]  *  2  )  +  (  6  xG:.'/[  3  ;  ]  X  ^C  2  ]  ) +  C;/C  4  ;  ]  xG:[  1  ]  *H 


Cl] 


7 
7 

7 
7 


Z-^y    VATAN    X;C 

Z-H((0>^)x(ol)xo  =  Y)t((  xy)xCxoi2  )  +  (~C)x(  (xj)x(o>;i:)xoi)  + 

~3oYiXi-C'^X=0 


[1]  Z^l 

[2]  -^(l=p,S)/8 

[3]  ->((0    1    xpS)=/;^"l  +  (0=S[  ;l])i0)/    6    9 

[4]  /JCl+A',0;]-<-Z?ri  +  o,,V;] 

C5]  SCl;  ]-H-Bni;] 

[6]  B-t-Cl    1    +S)-(1+S[  :l])o  .xl+£?[l  ;]i;/^l  +  ,B 

[7]  ->-2,Z-^Z^N 

[8]  -^OxpZ-«-Zx  ,5 

[9]  -►Z-^O 
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7  Z-^N   QD    n 

Cl]  Z^((2t//,0)  QP    Dill,-DZ22)i,liD    QP    D[i:i,-DL2l 
V 

V  Z^X   QP   Y 

[1]  Z^i-/Xxy),+/Xx(pY 
7 


V  Z*-N   QR    XiTiR 
[1]    2'*-(o2)|^C2]  ATAN    ,Y[  11 
[2]    B'^(+/X*2)*0.5 
[3]    Z^(;?*i//)xa)(2,//)p(20?'),lor-H(i;!/)xT+(o2)x~l+iiV 

7 

7    T   LIWA'3    L;KiUiAi3iCiD:XiYiZ;HinDiXDiyDiZDiXDD;yDDiZDD 
lZZ;DD',ND 

[1]       ;l'^(p2')pO  ' 

[2]         ;lCA';2]-e"lo(loL[l])xioi'[A'-^-    1    4    ;1] 

[3]         /l[^;l]-^((10L[l])x(202'[;/;l])f20/lC//;2])    VATAfK 

20LLll)i20AlNi2l 
[4]         /lC2;2]-^r[2;l]x(10L[l])x(202'ri;l])T20/l[l  ;2] 
Cs:         >l[2;l]-K-rC2;l]x(3oL[l])x(io2'[l;l])x(2o;l[l;  1])* 

2 
[6]         /l[3;2l"H(((loLCl])x(rC3;l3x2oT[l;l])-x/rC2;    1    m.lorCl 

;!])  +  ( lO/5ri;2])xyl[2;2]*2)^2  0>lCl;2] 
[7]         /l[3;l]^(3oL[l])x(x/2.2'C2;l],/l[2;l],    1    2    1    orC  1  ;1  ]  ,/l[  1 

;     1     l])-(2x.O/lCl  ;     1     l])x(r[3;l]xlor[l;l3)  +  x/J'[ 

2;    1    l],202'[i;l] 
[8]         /1[1    1    ;3]-h2'[1    4     ;2]+i?*-((20L[l])xlorf  1    U     ;1])    VATAf! 

2o!rCl    4    ;1] 
C9]         >l[2;3]-^rC2;2]+i?D-.-!rC2;llxB-H(20L[l])x((20/?^lti?)  + 

2or[l;l])*2 
C 1 0  ]      /I  [  3  ;  3  l-^rC  3  ;  2  ]  +  (  rr  3  ;  1  ]  xy,'  )  +  2  X  fC  2  ;  1  ]  X  (  20LC  1  ]  )  X  (  ♦  / 

20i?,7'Cl;l]  )x((r[2;l]x(  20i?)xi02'[l;l])-i?Z)x(10.'7)x 

2orCl;l])*(2orCl;l])*2 
Cll]    LO'.CHiF-^-ZZ-^-Ailill    AB    /I  C  1  ;  2  ]  ;  C/? 
Cl2]      i?-»-0    120    IfSP   A 
[13]       /?-H.'?,DF(;x(~20+/.'?[  ;i3]x,?[  ;     4     5    6  ]  )  .  (  ~202Z+ .  x<5/?[  ;  l 

3]). [1.1]    "2oZZ+.xa}/f[  ;     4    5     6] 
[14]        3pCn 
[15]       -►0xi2*p,,V 
[16]      ;f-^(loL[2])xlor[-V;2] 

[17]  y^-(l    2    x.oI)  +  (2  1  x.oi)x2or[/'/;2] 
[18]   Z-^(2x.OL)-(lx  .OZ,  )x20f[.V;2] 
[19]   4[.'/;2]-H~lo(Jx2o2'[//;l])-yxio2'[A';l] 
[20]   i4[/7;l]-^((-(Ixior[//;l])i-yx202'[//;l])  J20>1[A';2])  Ki^Tyl.V  Z* 

20/l[,V;2] 

[21]     zz?-^(y-t-i+y)  ,(z-«-i+;^)  ,z-t-itz 

[2  2]   .YD-«-3'C2;  2]x(loLr.  2]  )x2or[l  ;2] 
[23]   yP*-r[2;2]x(2  1  x  .oL)  xior[  1  ;  2  ] 
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[24]       ZD'*-Tl2;2l^(.l^  .oL)^loTlli2l 

[25]      A[.2i2l^(((2oniill)^XD-y^Tl2il^)-(loT[.l;ll)^yD  +  Xx?[ 

[26]       //^((202'Cl;li)xn^(yxZ/?)-Zxjn+;fxT[2;l])  +  (102'[l;l])xC^(Z 

xZD)-ZxXD-y^Tl2lll 
[27]       /l[2;l]-^^/i0^,'!x2T(2O/3ri  ;1])*2 
[2  8]       A:DZ?^(rr3;2]x(loLC2])x2or[l;2])-(7'[2;2]*2)x(loL[ 

2])xlor[l;2]  ..     ,     r       .-, 

[29]       YDD^{2    1    x.oL)xZDP^(7'C3;2]xlorCl;2])  +  (2'C2;2]* 

2)x2o2'[l;2] 
[30]       ZDDiril^.oL)^0  +  7.nD 
[31]      /l[3;2]-((2or[l;l])xIi:)/?-(yx2'[3:l])  +  (Ixr[2;l]* 

2)  +  2xyDxr[2;l])-(lo7'ri;l])xy/^-'^+(2xJZ)xr[2;l])  +  (XxT[ 

3;l])-yxr[2;l]*2 
[32]       ;l  [  3  ;  2  X  4  r  3  ;  2  ]  +  ( /I  C  2  ;  2  ]  *  2  )  X 1  o  /.  [  1  ;  2  ]  )  V  2  o  /I  [  1  ;  2  ] 
[33]       WD^(2or[l;l3)x(Cn'}xj)-(ZnxA'x7[2;l])  +  (-x-rPxr[ 

2  ;  1  ]  )  + (  2xXx5:>[  3  ;  n  ) +  (  Zx  y^5  ) -'^xyi:  2  ;  1  ] 
[34]       ,"/£)^/;D  +  (102'ri;l])x(ZDDx;f)  +  (Zi;xyxrr2;l])+(ZxYZ?x2'[ 

2;l])  +  (2xyxr[3;l])-(ZxV9p)+,'?xrC2;l] 
[3  5]      DD^(2xZi(20A[l;l])*2)xZi'^+4[2;l]xZx30/l[l;l] 
[35]      /[3;l]^((IJxA'.'^)-0Dx;/)^Z;x,r)  ,     ,,    , 

[37]       P-HC/Z-^dx.oDxlorCl    4    ;2])    y/ir/l.V    <':-»-(2oL[1])-( 

20L[2])xx/204ri    4     ;     1    2] 
[38]       All    n    i3-]^TLl    4     ;3]  +  P+.'7^((104[l    4     ;2])x20/l[l 

4     ;!])     7^7';4/7    lO/Ul    U    ;l] 
[39]       A?D^(lx.0Z,)xr[2;2]x20l'Cl;2] 
[40]      DD^(2oL[2])x(;l[2;l]x    i    2    x.o4[l;    1    2])+>l[2;2]x 

2  1    x.o/lCl;    1    2] 

[41]       PD«f-((2oP-Hl+P)*2)x((Z)x;//})-Z?Z?x//^l1-/V)TDx/)^ltZ? 
[42]      /l[2;3]--!r'[2;3]  +  PP  +  5P-((205^1+5)*2)x(/ir2;2]x    2 

3  *.o4Cl;    2    l])-/l[2;l]x(lo^Cl:2])  +  (lo/lCl;l])* 

2 
[43]  NPD^i  1  X  .  oL)  X  (  ?[  3  ;  2  1 X  20I'C  1  ;  2  ]  )  -  (  2'[  2  ;  2  ]  *  2  )  X  lore  1  ; 

2] 
[44]       Z)ZP/7*-(20ir2])x(/l[3;l]x     1    2    x.o^fl;    1    2l)i-(A[.3; 

2]x    2    1    x.o^Cl;    1    2])  +  ((.4C2;    1    2  ]  +  .  *2  )  x  2  x  .  o4  [  1  ; 

1  2])-'</2,4[2;     1    2l,loAll:     1    2] 

[45]       P-  (  -  2x  PDx  PD  X  jop  )  +  ((20P)*2)x((Dx  (.D^tJOP  )  -N^DDD  )  - 

2'^0I)x(Px»;") -A'x/7.0)  vD*3 
[46]       ^t3;3]^r[3:3j+?t(-2x5DxJ/')x305)  +  ((205)*2)x((20/l[l; 

2  ]  )  X  (/![  3  ;  2  ]  T  30/.  [  1  ;  1  ]  )  -  2x4[  2  ;  1  ]  xyl  [  2  ;  2  ]  t  ( lo/l  [  1  ;  1  ]  )  * 

2) 

[4  7]  4[3;3]^/l[3;3]  +  ((20j)*2)x(104ri;2])x(2x(/l[2;l]* 
2)t  1  1  3  x.O/5[1;1])-(/1C3;1]t(10/,[1;1])*2)  +  (4[ 
2;2]*2)i304[l ;1] 

[48]       -^LO 
.V 
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